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INTRODUCTION 


The intent of this manual is three-fold. It is, first of all, in- 
tended to be an introduction to COBSTRAN (Composite Blade STRuctural 
ANalyzer) by providing an overall description of the program capabilities 
and its interaction with general purpose structural analysis codes. 
Secondly, it describes the organization of the program and the flow of 
data and analysis through the subroutines. Thirdly, it describes the 
equations formulated and solved, and the function of each of the 
seventy-two subroutines. However, the detailed description of the theory 
used in developing the programmed equations, especially in the composite 
mechanics and laminate analysis procedures, is properly left to the ref- 
erences from which they were obtained.. 

Other manuals available are a USER'S MANUAL, describing in detail 
the input deck set-up and options available to the user, and a DEMON- 
STRATION MANUAL with selected problems. 

COBSTRAN is designed to carry out the many linear analyses required 
to efficiently model and analyze blade-like structural components made 
of multilayered angle-plied fiber composites (Figure 1). 

Inputs to COBSTRAN are composite constituent fiber and matrix mate- 
rial properties, factors reflecting the fabrication process, composite 
geometry and model geometry. COBSTRAN determines the ply layup at each 
grid point and calculates the equivalent homogeneous properties at each 
grid point. 

Outputs are the individual ply properties and nodal composite prop- 
erties consisting of membrane, bending and coupling stiffnesses. Also, 
finite element structural analysis grid, connectivity and material data 
are generated. COBSTRAN generates model finite element data in formats 
compatible with COSMIC NASTRAN, MSC/NASTRAN and MHOST. MHOST is a Lewis 
Research Center structural analysis code based on mixed-iterative finite 
element methodology and is specifically designed for three-dimensional 
inelastic analyses. Reference to the MHOST code in this manual is for 
completeness only and its use is discouraged until release of the MHOST 
code and a fully tested COBSTRAN version. 

The laminate decomposition analysis procedure is designed to make 
use of element stresses from the structural analysis code and information 
generated in the COBSTRAN preprocessor phase which is stored in temporary 
files. Membrane forces and moments are calculated in the structural co- 
ordinate system for each node in the discrete model of the structure. 

For each ply layer at each node the stress and strain state is calculated 
including the results of combined-stress strength failure criteria and 
Hoffman's failure criteria (Ref. 1). 
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An additional feature of COBSTRAN Is that the geometry input may be 
an existing COSMIC NASTRAN or MSC/NASTRAN bulk data deck. Material prop- 
erties are generated within COBSTRAN and a complete finite element model 
is created in any of the optional output formats. For more details see 
the NASMIN and MSCMIN program option cards in Section 3 of the USER'S 
MANUAL. 


1-2 (8-31-88) 



COBSTRAN ITieoretical/Programmers Manual 


COMPONENT 


FINITE ELEMENT 
MODEL 


LAMINATE 



♦ 


LAMINATE 

THEORY 



CONSTITUENTS 



LAMINATE 

THEORY 


♦ 




PLY PROPERTIES 
a e F.C. 


Figure 1 . - COBSTRAN analytical process. 
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SOLUTION FORMATS 


Two types of problem solution formats are available in this code. 
The first format is for a solid composite blade or solid composite blade 
with hollow core sections of the turboprop or turbofan types. Plate or 
shell type 2-D elements are created along the blade mid-thickness line 
and equivalent homogeneous anisotropic material properties are calculated 
from the composite layup for each element. The second format is for 
blades of the aircraft wing type construction composed of quadrilateral 
shell elements on the airfoil surface separated by plate type spars at 
user selected chord positions. This format is most suitable for long 
wing-like structures used for wind power generators. This code would also 
be suitable for the analysis of any composite structure which can be 
modeled using plate bending or membrane type elements. 

As a preprocessor only, COBSTRAN can be used to generate 
nodes(grids), connectivities and material properties for various config- 
urations which can then be combined into one structural analysis model. 
COBSTRAN writes these bulk data records to FORTRAN unit KBULK which can 
be assigned to a user designated file. COBSTRAN was designed primarily 
to interface with the NASTRAN structural analysis code, therefore, some 
of the terminology used in this manual will be familiar to NASTRAN users. 

When assembled the COBSTRAN input data deck consists of the following 
five parts. 

1. Option deck 

2. Parameter deck 

3. Composite deck 

4. Geometry deck 

5. Composite databank (in separate file) 

The Option deck identifies the job and type of solution to be per- 
formed. Twenty user-selected options are available establishing such 
items as model input/output format, composite ply order and preparation 
of data input for the general purpose structural analysis code to be used. 

The Parameter deck defines geometry and material property 
input/output characteristics such as mesh density, number of composite 
material systems (i.e. fiber/matrix combinations), ply stacking sequence 
and type of material output (membrane, bending, etc. ) for use by the 
structural analysis code. 

The Composite deck defines the characteristics of each ply layer. 
These are fiber type, matrix type, thickness, fiber volume ratio, orien- 
tation angle and ply location, by percent chord and percent span, in the 
structural model. 


1-4 (8-31-88) 



COBSTRAN Theoretical/ Programmers Manual 


The Geometry deck defines the model configuration with a minimum of 
input values. Four different input formats are available for this data. 
COBSTRAN 1 s internal interpolation routine is used to generate and refine 
the mesh pattern as specified by parameter values. 

An internal databank is provided with the COBSTRAN code which con- 
tains constitutive properties of 21 fiber types, 17 matrix types and 
correlation coefficients for 31 fiber/matrix combinations. This property 
databank is provided for user information and guidance only and no re- 
sponsibility is assumed by the authors as to accuracy or currency of the 
data. All properties are in customary U.S. units. The databank is stored 
in subroutines BANKCO, BANKFP, and BANKMP. 

For a user supplied databank, a format and definition of properties 
is provided in Section 4 of the USER'S MANUAL. This databank should be 
maintained in a separate file which is referenced to FORTRAN unit KDBANK 
when used with COBSTRAN. 

COBSTRAN provides a hierarchial approach to databank aquisition. 

a) The user supplied databank is searched for the 
fiber and matrix properties requested. 

b) The internal databank is searched for those 
properties not found in the user supplied 
databank. 

The source of the properties is indicated in the output. 


MESH GENERATION 


The mesh generation capability in COBSTRAN is designed to create 2-D 
plate type meshes under the SOLID option and aircraft wing type meshes 
under the SHELL/SPAR option. In both cases meshes are created from 2-D 
elements. Under the SOLID option four types of 2-D mesh patterns may be 
created from nodal points generated along the mid-thickness line. These 
patterns are described in User Controlled Parameters, Section 3 of the 
USER'S MANUAL. 

Under the SHELL/SPAR option only quadrilateral elements are uti- 
lized. Each airfoil chordwise cross-section is divided into approximately 
equal spaces but modified for grid points at spar locations. The leading 
edge section of each airfoil surface is divided into two elements for an 
improved mesh on a small radius surface. 
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COMPOSITE MECHANICS 


The properties of each layer at each node are calculated using the 
micromechanics equations programmed in subroutine PLY. These individual 
layer properties are combined using the laminate theory equations pro- 
grammed in subroutine PROP forming the elastic property matrix for the 
node. PLY and PROP are cycled for each node. 


LAMINATE DECOMPOSITION 


The element stress or strain output from the structural analysis code 
is used in subroutine STR to calculate the two-dimensional load vector 
at each node of the descrete structural model. Using these load vectors, 
the subroutines PLY, PROP and COMSA are cycled for each node calculating 
the stress and strain state and failure criteria for each lamina. 


OUTPUT DATA 


The laminate analysis module generates material properties for each 
node. The nodal values printed are; 

Membrane stiffness matrix 
Bending stiffness matrix 
Coupling stiffness matrix 


The finite element data module generates and writes to FORTRAN unit 
KBULK the following; 

Element connectivities 
Nodal coordinates 
Element material properties 


The laminate decomposition module generates and reports for each 
node the stress, strain and failure criteria for each ply at the node. 
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A summary of the minimum margins of safety for each node is printed 
as follows; 

Combined stress margin of safety 
Hoffman failure margin of safety 
Interlaminar shear margin of safety 


1-7 ( 8 - 31 - 88 ) 






COBSTRAN ITieoretical/Programmers Manual 


INTRODUCTION to Program Modules 


This section describes the function of each subroutine in COBSTRAN 
and where appropriate, the equations programmed in each subroutine. The 
flow of COBSTRAN can be followed by examining the flow chart of the MAINPG 
program. There are four secondary level executive subroutines called in 
MAINPG. These are CSOLID to generate a two-dimensional model, CSHELL to 
generate an aircraft wing type model, SOLVEQ to drive the MHOST program 
and CPOST to analyze the output of each finite element analysis program 
and calculate the individual ply stresses, strains and failure criteria. 

A list of engineering symbols with their FORTRAN code and description 
are in SYMBOLS and Code Identifiers, Section 4 of this manual. Additional 
FORTRAN variables and their descriptions are available in the DICTIONARY, 
Section 6 of the USER'S MANUAL. 
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MAINPG 


CALLS TO — CINIT 

REFLEC 

LOGO 

OPTION 

CSOLID 

CSHELL 

SOLVEQ 

CPOST 


COMMON 

/CUNITS/ KDES1, 

KKEI, 

KBULK, 

KREAD, 

KOUT, 


KTTAB, 

KDES2, 

KPROPI, 

KGNP, 

KDBANK, 


KPR0P2, 

KPLOT 

KPOST, 

KMHFOR, 

KREFL, 

KDIAG, 

COMMON 

/BOOL/ BIDE, 

DYNA, 

PRTOUT, 

WINDML, 

SOLID, 


MINICH, 

UNSY, 

PROFIL, 

PROT, 

PREP, 


POST, 

MODES, 

HOLLOW, 

COSD, 

MSCD, 


PLYORD, 

RINDV, 

TEM, 

ZIGZAG, 

PRES, 


COSMIC, 

RANG, 

MSCNAS, 

DUMP, 

MHOST, 

CONE, 

EZRD, 

MPLOT 

ECHO, 


COMMON /DIAGS/ LDIAG(IOO) 


The main program initializes program FORTRAN units, stores the us- 
er's input deck for later reading, generates the logo and sets program 
options. The options selected are used to control program flow through 
the four second level executive subroutines, CSOLID, CSHELL, SOLVEQ and 
CPOST. The flow charts presented in CSOLID, CSHELL and CPOST describe 
the subroutines utilized within each of these paths in COBSTRAN. 

To modify this program for use with the MHOST code remove 1 CHOST 1 
from the lines so labeled in this main program. Subroutine SOLVEQ at the 
end of this code is the main program for the MHOST analysis code. 

These codes are setup for large word machines -64 bits-, to change 
to double precision, remove 1 CDP 1 from the appropriate lines in this code 
and modify subroutine SOLVEQ according to its instructions to change the 
MHOST code. 
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SUBROUTINE BANKCO 


CALLED FROM — BANKRD 


CALLS TO — none 


ARGUMENTS (NA1,NA2,C0R, IF, IX) 


Correlation coefficients for the following fiber/matrix combinations 
are stored in this subroutine. 


Fiber/Matrix Date Comment 


T300/SPOX 

T300/EPOX 

ALT6/ALT6 

NIKL/NIKL 

HGLA/EPOX 

FOAM/FOAM 

HT-S/GK60 

MOD1/EPOX 

M0D2/EP0X 

M0D2/P0LY 

HM-S/EPOX 

HT-S/EPOX 

MODI/POLY 

HT-S/POLY 

B0R4/EP0X 

SGLA/EPOX 

KEVL/EPOC 

HV-S/GV60 

HT-S/PMRP 

B0R4/ALT6 

KEVL/TIT6 

MOD2/BERY 

BOR5/EPOX 

KEVL/EPOX 

SW4M/NM40 

HERC/SSAL 

BOR5/ALT6 

TITF/TIT6 

ADHF/FMTH 

HERC/EPOX 

HOLO/HOLO 


9-26-80 FOR COMPOSITE PROPERTY INPUT BY USER 

10-26-79 
8-01-83 
8-12-83 

10- 11-83 CROSS WEAVE FOR SR7 PROP MODEL 
8-08-83 

11- 27-72 
06-15-72 
06-15-72 
06-15-72 
06-15-72 
06-15-72 
06-15-72 
06-15-72 
11-07-73 
11-07-73 
11-12-73 
11-27-72 

06- 20-74 
11-07-73 

07- 10-74 
07-10-74 
02-04-75 
02-04-75 
02-10-75 

02- 04-75 
01-25-77 
01-25-77 
01-25-77 
01-25-77 

03- 20-79 
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SUBROUTINE BANKFP 


CALLED FROM — BANKRD 


CALLS TO — none 


ARGUMENTS (NAF,I1,F1,I2,FP,IF,IX) 


Properties of the 
routine. 


Fiber 

Date 

ALT6 

ISO 

8-01-83 

NIKL 

ISO 

8-12-83 

HGLA 

ISO 

10-11-83 

FOAM 


8-08-83 

T300 

ANIS 

10-26-79 

M0D2 

ANIS 

06-15-72 

MODI 

ANIS 

06-15-72 

HM-S 

ANIS 

06-15-72 

HT-S 

ANIS 

06-15-72 

AS— 

ANIS 

06-15-72 

B0R4 

ISO 

11-07-73 

SGLA 

ISO 

11-07-73 

KEVL 

ANIS 

11-07-73 

HV-S 

ANIS 

06-15-72 

BOR5 

ISO 

04-11-74 

B0R8 

ISO 

04-11-74 

EGLA 

ISO 

04-11-74 

SW4M 

ISO 

02-10-75 

TITF 

ISO 

01-25-77 

ADHF 

ISO 

01-25-77 

HOLO 


03-20-79 


following fiber materials 


Comment 

NO REF 
NO REF 

CROSS WEAVE FOR SR7 PROP 
FOR SR7 PROP MODEL (SR7A 
NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 

NO REF. 


are stored in this sub- 


5-23-84) 
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SUBROUTINE BAN KM P 


CALLED FROM — BANKRD 


CALLS TO — none 


ARGUMENTS (NAM, II , PM, IF, IX) 


Properties of the following matrix materials are stored in this 
subroutine. 


Matrix 

Date 

SPOX 

ISO 

9-26-80 

GK60 


11-27-72 

FOAM 


8-08-83 

NIKL 


8-12-83 

EPOX 

ISO 

06-15-72 

ERLA 

4617 

06-15-72 

POLY 


06-15-72 

ALT6 


11-07-73 

TIT6 


11-07-73 

EPOC 


11-12-73 

SSAL 


11-12-73 

GV60 


11-27-72 

PMRP 


06-20-74 

BERY 


07-10-74 

NM40 


02-10-75 

FMTH 


01-25-77 

HOLO 


03-20-79 


Comment 

FOR COMPOSITE PROPERTY INPUT BY USER 
(ORIGINAL MODEL SR2C BLADE) 

FOR SR7 PROP MODEL (SR7A 5-23-84) 

NO REF. 

NO REF. 

NO REF. 

MAT. ENG. 1967 
NO REF. 

MAT. ENG. 1967 


MFCA DATA 
NO REF. 
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SUBROUTINE BANKRD 

CALLED FROM — CSHELL CSOLID 

CALLS TO — BANKCO BANKFP BANKMP 

ARGUMENTS ( CODES, FP,MP,CORC,DF,NF, NT, NPD,NDE$,KDBANK,KOUT) 


Composite material system descriptions are read into the program by 
this subroutine. Fiber properties, matrix properties and correlation co- 
efficients from a user supplied databank must be assigned to each COBSTRAN 
preprocessor run. This databank is stored in FORTRAN unit KDBANK. 

Subroutine BANKRD first searches the user supplied databank for the 
fiber and matrix properties selected in the input deck. If properties 
are not found in the user supplied databank the internal databank residing 
in subroutines BANKCO, BANKFP, and BANKMP is searched for the missing 
properties. Therefore, if only the properties from the internal databank 
are to be used, a user supplied databank is required containing three 
records with zero in the first column of each record. Each record refers 
respectively to fiber, matrix and correlation coefficients. 

An output message informs the user regarding the source of the 
properties. 

The selected fiber and matrix properties and the corresponding cor- 
relation coefficients are retained in a temporary file for 
postprocessing. 
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SUBROUTINE CINIT 


CALLED FROM — MAINPG 
CALLS TO — none 


ARGUMENTS — none 


COMMON 

/CUNITS/ KDES1, 

KKEI, 

KBULK, KREAD, 

KOUT, 



KTTAB, 

KDES2, 

KPROP1, KGNP, 

KDBANK, 



KPR0P2, 

KPLOT 

KPOST, 

KMHFOR, KREFL, 

KDIAG, 

COMMON /BOOL/ 

BIDE, 

DYNA, 

PRTOUT, WINDML, 

SOLID, 



MINICH, 

UNSY, 

PROFIL, PROT, 

PREP, 



POST, 

MODES, 

HOLLOW, COSD, 

MSCD, 



PLYORD, 

RINDV, 

TEM, ZIGZAG, 

PRES, 



COSMIC, 

MSCNAS , 

MHOST, EZRD, 

ECHO, 



RANG, 

DUMP, 

CONE, MPLOT 


COMMON 

/WINDS/ 

SY(4,30,5) 

, MWSCT(30,2), ALPHA(3,3) , 

X(30), 



NX$PC(30), 

NSPDES(5), NGRID, 

NELT, 



NPAST, 

THMAX, 

NXSECT, 

NSPAR, 



NXSPAR, 

NYSPC, 

LSP(2,5,30), 

NYPTS, 



NNQ, 

NNX, 

LSECT(2) , 

NPTOT, 



NETOT, 

HOLD(IO), MAXPLY, 

MPLY(300) 


This subroutine establishes the FORTRAN units to be used in COBSTRAN. 
If the MHOST option is selected a second call is made to this subroutine 
to reset the units to avoid conflict within the MHOST module. 


Variable 


FORTRAN 

FORTRAN unit 

name 


unit 

for MHOST 

* KDES1 

- 

52 

52 

* KKEI 

= 

53 

53 

* KBULK 


54 

54 

KDIAG 

= 

6 

6 

KREAD 

= 

5 

5 

KOUT 

= 

6 

6 

KTTAB 

= 

57 

57 

KDES2 

= 

58 

58 

KPROP1 

r: 

59 

59 
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KGNP = 

60 

60 

KDBANK = 

61 

61 

KPR0P2 = 

62 

62 

KREFL = 

63 

63 

KPLOT = 

64 

64 

KPOST = 

5 

55 

KMHFOR = 

90 

12 


* Data stored on units KDES1, KKEI, KBULK and KGNP must be saved for 
post-processing of the general purpose structural analysis code. All 
other units are scratch files. 
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SUBROUTINE COBSMF 


CALLED FROM — COMSA 
CALLS TO — none 

ARGUMENTS ( J , FPT, MPT, SIGL, KSMF, SMF, PL) 


The strain magnification factors from which the ply unidirectional 
limiting are constructed are generated in this subroutine. These factors 
are 4 > 22 * ^12 anc * ^23 » f° r constructing anc * S^23 respec- 

tively. 

Four user-optional methods are available to compute 0^2 » Parabolic 
strain distribution, Daniel’s indirect, Kies's two-dimensional , and Linear 
strain distribution. 

Filament and matrix orthotropy and the effects of voids are included 
in all of these methods as described in reference 2. 

The equations programmed in this subroutine are; 


VC 


1.0 


4k v \ 1/2 
* k m ) 


k f = (l-k v ) k f 


k m ~ ^ k f) 


k m — (1 ~ ky)(l k-p) 


C ^ 


^m^ ~ 


(1 - kf) 
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The four methods used to compute the strain magnification factor 
</> 2 2 » are as follows; 


(1) Parabolic strain distribution method; 



Em22 

E f22 


2 PR 

+ PR- 1.0 


+ 


1.0 


J 


**22- PL(43,J) 


(2) Daniel's indirect method 


«V22 = M 1 - ^23) 


e ^12 

E ml2 


O. 83 IY -^- V /2 -2 

2 - 1 . 3s[ f 

” ) 

|^ 2 - 2 

[\kf J \ 


k f ) 


0.35< k f <0.75 

if not 

then 

**22 = 


+ 1.78 


= 1.0 


V2( VC >- PL(44,J) 
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(3) Kies's two-dimensional method 


_ ^ (1 - v f 12 v f21 ) E m22 
t 1 - v nl2 v ro2l) c f<#> E f22 

E = , ml2* 


SI = -f— - C 1 - V21 P <'fl2 “ B >] "(ZZ 
V 22 


S2 = -f— [P (»f 12 - B> - VI 2] "el\ 
Vll 


V 2 


m22 = 

1 

1 + p (A 

-1) 

-j(Sl + S2) 


_ V22 


•Vll 

V22 

e ^22 

VI 2 C „ 
*V11 

£ m22 


if 1 

V22 1 a l c 

V22 



1.0 


if- 

£m22 <1.0 
V22 “ 


,-30 


*„22< vc >- PL(45, J) 
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(4) Linear strain distribution method 


^,22 


PR 

PR- 1.0 

E f22 


*,22 (VC)- PL(46,J) 


STRAIN MAGNIFICATION FACTOR ^ 12 


*«12 = 


1-P 1- 


E rol2 
c fV E f 12 


*,,12 (VC)-. PL(47,J) 


STRAIN MAGNIFICATION FACTOR </. 2 3 


23 


2 (1 - P) + (2p -1) 


E m<(> E m23 
c f 4 , E f23 


O 3 (VO- PL(48,J) 
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Inputs to subroutine COBSMF are the ply applied stresses (c^n, 

V22 * V33^» the volc * anc * apparent fiber content, and ply fiber and 
matrix elastic constants. 

The outputs of this subroutine are the magnification factors stored 
in array PL(43,J) thru PL(48,J). The strain magnification factor selected 
by the user is transferred by the call argument SMF. 
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SUBROUTINE COMSA 


CALLED FROM — CPOST 

CALLS TO — COBSMF INVA PLYIN 


ARGUMENTS ( I , PROT, PRTOUT, RINDV, SPARS , KTTAB, KOUT, KSMF, 
MESH , NPAST , NLEP , FC , NSB , MSB , 

CODES , FP , MP , CORC , NT , NPD , THERMO , MDES ) 

COMMON /PROPER/ PL(71,223) 

COMMON /STRA/ DISV1(6,510) 

COMMON /PRINT1/ ID(80,5), IOPT(80,30) 

COMMON /PRINT2/ NPRT(2,10) 

COMMON /PLIN/ FPT( 19) ,MPT( 15) , CC(6, 5) ,TS3(3) , ALPP(3) 

COMMON /PROPTY/ EP( 16 ) , THP( 5) , THF( 6) , AS(9) , CS( 9 ) , 

BS(9) ,RA(9) , RB(9) ,DENS 


In this subroutine the stress and strain state of each ply are com- 
puted given the edge membrane force resultants, the ply temperature change 
and the changes in curvature, In addition, two-ply combined-stress 
strength criteria and the interply delamination criterion are generated. 
The equations programmed for the i^h ply strain and stress states are 


{Vi> ~ CR/i - 1 [A CX ] + {N cA j x } + ^cx-H w cbxO zCR^*] {w c b x } 


(Vii-CE/i] *£fyl^EA cx ] ^ ({N cx } + {N c aTx^ + ^cx^ w cbxO 

“ ^'i ^ + z ^/i ^ ( w cbxO 
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The composite reference plane strains c csx and the changes in cur- 
vatures are computed from 


j { £ csx} 

[ { w cbx* 


[A CX J | CC cx 3 
[C cx ] I tD cx ] 


-1 


f ^cx^ I f < n caTx> V ] 

( ^cx^ J [ ( m caTx) Jj 


when one or both of the membrane force and moment resultants are given. 

The strains are generated locally in EPSL and SIGL, respectively, 
and are stored in PL( 64, 1 ) to PL(69,I). The matrices CR n *] and [E^O are 
generated locally from information transferred from PL( 14 , I ) and PL( 3 1 , I ) 
to PL(42, 1). The distance z n - and the ply temperature change aT / 1 - are 
transferred from PL(11,I) and PL(50,I), respectively. 

The remaining matrices are; 


A cx - ACX 
^cx 

n caTx - NSDT 
N cx - NSB m 

M C x - MSB m 


w cbx WXX m (local curvature from bending analysis) 
where m denotes the load direction. 

It is important to note that the stress analysis in the coded form 
also handles the case where both reference plane membrane strains and the 
local curvatures are given. In this case the ply strains are given by 


(Vxi} ~ { £ csx} “ z e\ { w cbx^ 


where { £ ^ xl -} are the ^ ply strains along the structural axis, {f; csx } 
are the composite reference plane membrane strains, z^-j is the distance 
from the reference plane to the centroid of the ith ply, and are 
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the composite local curvatures. They are read in the array D vm where m 
denotes the load direction. 

The corresponding i^ ply stresses are given by 


{o-^i) - [E^i] ([R^i]{£^xi) A ^Vi 

where {o^} are the i^h ply stresses along the material axes, [E^] are 
the i^h p]y elastic constants about the material axes, [R^] is the 
transformation matrix of the i^h ply from the ply material coordinate 
system to the structural coordinate system, {^ X1 *} are the i^h ply 
strains along the structural axes as given by a previous equation, aT^ 
is the temperature change of the i^h ply, and {«,.,■} are the thermal co- 
efficients of expansion of the i^h ply i n its material coordinate system. 


The displacement force relations are shown; 


Di spl acement 

Di spl acement 

force relations 

Forces 

| JU C x> 1 


i — i 
X 

u 

< 

1 1 

tc cx ] ] 

-1 

f ^cx^ 

1 <«cx> J 


[C cx ] 1 

Wcx] J 


[ ^cx^ 


Output is by node(grid point) presenting ply stresses, strains and 
failure criteria for each ply layer. 
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The failure criterion may be determined by either of the following 
methods. 


(1) Modified distortion energy 


F = 1 - 


Vila \2 

*tiu ) 


a tzza \2 K vn« °az 
hzz? ) ~ ' W 'W 


V12S \2 . 

S/12S ) 



P L( 62 , 1) 


The parameters « and p are specified as follows: 


T Vll^° 

C vil <0 


{ T V22 ^ 0 

C v 22 < 0 


s ALla = 


S jfllT 
ml n(S^»nc» 


S ^11CD) 


a = T 
a = C 


K nZa(i = K Vl Zap 


(1 + 4v^ 12 - V13)E,22 + i 1 - V23^ E /11 
[ E ^ll e /22 i 2 + VI 2 + V13> (2 + V21 + V23>] 


1/2 


2-18 ( 8 - 31 - 88 ) 



COBSTRAN Theoretical/ Programmers Manual 


K Vl2 cep - 


/*C0RC(1,5) 
J C0RC(2,5) 
) C0RC(4,5) 
\.CORC(3,5) 


a, /? = T 
a = C , /? = T 
a = T, /?=C 
a, p = C 


The multi plyer of K^a/? was generated in subroutine COMSA and is 
stored in PL(61,I). The constants K ^ 12aJ? constitute theory-experiment 
correlation factors. 


(2) Hoffman's criterion (ref 3): 


VllC *" ml n ( s ^llC> s mCD> 


F= 1 - 


2 2 c 

VI 1~ VII V22 + V22 + V11C- VllT 

S ^11C 1 IT S «?22C S ^22T S /11C S ^11T ^ 


F 

F 

F 


2 

s ^22C ~ s ,f22T , VI 2 
S <P22C s ^22T f S| 12S 


> 0 no failure 
= 0 incipient failure 

< 0 failure 



PL(71, I) 


2 - 19 ( 8 - 31 - 88 ) 



(X)BSTRAN Theoretical/ Programmers Manual 


The interply delamination criterion for the jth interply layer at 
the m^h load condition is governed by 



/ I A</> I \ 

V A ^del r 


PL(63, I) 


when i > 1 


A</>j = -“ (£ C yy ~ £ CXX^ S1 * n ^ 0 i ~ S1 n ^^i-1^ + \ c CXy ( C0S ^i c °s20j_l) 

{£ CX } = tA cx l + {N cA j x } + [C cx l {w cbx }) 

or is given by the displacement force equation described previously. 


The inputs to the subroutine are the ply angle measured from the 
structural axes from PL(14,I), the distance from the reference plane to 
the centroid of the ply z n * from PL( 11 , I ) , the ply temperature change 
aT^ from PL( 50,1), the interply delamination limit A<t>d e -\ j from 
PL(60,I), and the ply thermoelastic properties stored in PL(24 to 26,1) 
and PL( 31 to 42,1). The ply extensional and coupling rigidities 
A cx = ACX and C cx = CPC; the local curvatures w cbx = WXX; the adjustment 
constants K'^2TT = C0RC(1,5), K'^j2CT = C0RC(2, 5), K'^2TC = C0RC(4, 5), and 
K Vl2CC = C0RC ( 3 > 5); and the load conditions N cx = N$B(m). 

The subroutine outputs are the modified distortion energy PL(62,I), 
Hoffman's criterion PL(71,I), the interply delamination criterion 
PL(63, 1) , and the adjacent ply relative rotation (A</>j from PL( 70 , 1 ) ) . 
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SUBROUTINE CPOST 


CALLED FROM — 
CALLS TO — 


MAINPG 

COMSA MHSTR 

PLYSHR PROP 


OUTPUT PLY 
STR 


ARGUMENTS none 


COMMON 

/CUNITS/ 

KDES1, 

KKEI, 

KBULK, KREAD, 

KOUT, 



KTTAB, 

KDES2, 

KPROP1 , KGNP, 

KDBANK, 



KPR0P2, 

KPLOT 

KPOST, 

KMHFOR, KREFL, 

KDIAG, 

COMMON 

/BOOL/ 

BIDE, 

DYNA, 

PRTOUT, WINDML, 

SOLID, 



MINICH, 

UNSY, 

PROFIL, PROT, 

PREP, 



POST, 

MODES, 

HOLLOW, COSD, 

MSCD, 



PLYORD, 

RINDV, 

TEM, ZIGZAG, 

PRES, 



COSMIC, 

MSCNAS, 

MHOST, EZRD, 

ECHO, 



RANG, 

DUMP, 

CONE, MPLOT 


COMMON 

/WINDS/ 

SY(4,30,5), MWSCT(30,2), ALPHA(3,3), 

X(30), 



NXSPC( 30) , 

N$PDES(5), NGRID, 

NELT, 



NPAST, 

THMAX , 

NXSECT, 

NSPAR, 



NXSPAR, 

NYSPC, 

LSP(2,5,30), 

NYPTS, 



NNQ, 

NNX, 

LSECT(2), 

NPTOT, 



NETOT, 

HOLD(IO), MAXPLY, 

MPLY(300) 


COMMON /VALUES/ NP, NE, NDES, IMAT, KSMF, IU, JU, XBEG, XEND, 


NMODE(3) , NOPLY, NSECT, MESH, NGROUP, NCOREG, 
IGRD, MSECT(40) , NC00R(3,8) 

COMMON /PLYS/ THERMO(223) ,MDES(223),THHF,NLEP 
COMMON /PROPTY/ EP(16),THP(5) ,THF(6),AS(9),CS(9), 

BS(9) , RA(9) , RB(9) ,DENS 
COMMON /PROPER/ PL(71,223) 

COMMON /MFC A/ C0DES(6,40), PERT(6,40), FP(19,1,40), 

MP( 15, 1,40) , CORC(6,5,40), DF(40), 

NF(40) , NT(2,40), NPD(40) 

COMMON /TABLES/ GNP(12,510),KEI(6,910),NCORDL(910),RAD(2,510) 
COMMON /STRA/ DISV1(6, 510) 

COMMON /PRINT1/ ID(80,5), IOPT(80, 30) 

COMMON /PRINT2/ NPRT(2,10) 


This subroutine is a second level executive program. It's purpose 
is to control the postprocessing phase of the analysis. Element stresses 
are read from the finite element analysis generated output file and 
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transformed to the COBSTRAN structural coordinate system. Nodal stresses 
and corresponding nodal membrane forces and moments are calculated. From 
these loads, ply stresses, strains and failure criteria are determined 
for each ply layer at each grid point. 


2-22 ( 8 - 31 - 88 ) 



COBSTRAN Theoretical/Programmers Manual 


DO LOOP 
I = 1, NP 


MHSTR 


1 S1 

rR 





PLY 



PROP 



COMSA 



-b- 




Figure 3. * CPOST flow chart. 


2-23 ( 8 - 31 - 88 ) 





COBSTRAN Theoretical/ Programmers Manual 


SUBROUTINE CSHELL 


CALLED FROM — MAINPG 


CALLS TO — BANKRD 

MHSHPR 

NASTDK NOLAPR 

OUTPUT 

OUTP11 

PLTBLK PLYIN 

WIND 

WINOUT 

WINPLY WINPUT 

WREAD 

WROTAT 

WSPAR 


ARGUMENTS none 






COMMON /CUNITS/ 

KDES1, 

KKEI, 

KBULK, 

KREAD, 

KOUT, 


KTTAB, 

KDES2, 

KPR0P1 , 

KGNP, 

KDBANK 


KPR0P2, 

KPLOT 

KPOST, 

KMHFOR, 

KREFL, 

KDIAG, 

COMMON /BOOL/ 

BIDE, 

DYNA, 

PRTOUT, 

WINDML, 

SOLID, 


MINICH, 

UNSY, 

PROFIL, 

PROT, 

PREP, 


POST, 

MODES, 

HOLLOW, 

COSO, 

MSCD, 


PLYORD, 

RINDV, 

TEM, 

ZIGZAG, 

PRES, 


COSMIC, 

MSCNAS, 

MHOST, 

EZRD, 

ECHO, 


RANG, 

DUMP, 

CONE, 

MPLOT 



COMMON /DIAGS/ LDIAG(IOO) 

COMMON /WINDS/ SY(4,30,5), MWSCT(30,2), ALPHA(3,3), X(30), 

NXSPC(30), NSPDES(5), NGRID, NELT, 

NPAST, THMAX, NXSECT, NSPAR, 

NXSPAR, NYSPC, LSP(2,5,30), NYPTS, 

NNQ, NNX, LSECT(2), NPTOT, 

NETOT, HOLD(IO), MAXPLY, MPLY(300) 

COMMON /VALUES/ NP, NE, NDES, IMAT, KSMF, IU, JU, XBEG, XENO, 
NMODE(3), NOPLY, NSECT, MESH, NGROUP, NCOREG, 
IGRD, MSECT(40), NCOOR(3,8) 

COMMON /PLYS/ THERM0(223),MDES(223),THHF,NLEP 
COMMON /PROPTY/ EP(16),THP(5),THF(6),AS(9),CS(9), 

BS(9) , RA(9) , RB(9) , DENS 
COMMON /PROPER/ PL(71,223) 

COMMON /MFC A/ CODES(6,40), PERT(6,40), FP( 19,1,40), 

MP(15,1,40), CORC(6,5,40), DF(40), 

NF( 40), NT( 2 , 40 ) , NPD(40) 

COMMON /TABLES/ GNP(12,510),KEI(6,910),NCORDL(910),RAD(2,510) 
COMMON /STIFF/ EX(3,510),AB(6,510),DEN(510) 

COMMON /STRA/ DISV1(6,510) 
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This subroutine is a second level executive program. It's purpose 
is to control the generation of aircraft wing type models composed of 
quadri laterial shell elements on an airfoil surface separated by plate 
type spars. 
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SUBROUTINE CSOLID 


CALLED FROM — MAINPG 

CALLS TO — BANKRD 

CSREAD 

MESHG 

MHPREP 


MSCMIN 

NASTDK 

NASMIN 

NOLAPR 


OUTPUT 

PLYIN 

0UTP10 

TBGRD 

0UTP11 

TBPLY 

PLTBLK 


ARGUMENTS none 

COMMON /CUNITS/ KDES1, 

KKEI , 

KBULK, 

KREAD, 

KOUT, 

KTTAB, 

KDES2, 

KPR0P1, 

KGNP, 

KDBANK, 

KPR0P2, 

KPOST, 

KMHFOR, 

KREFL, 

KDIAG, 

KPLOT 

COMMON /BOOL/ BIDE, 

DYNA, 

PRTOUT, 

WINDML, 

SOLID, 

MINICH, 

UNSY, 

PROFIL, 

PROT, 

PREP, 

POST, 

MODES, 

HOLLOW, 

COSD, 

MSCD, 

PLYORD, 

RINDV, 

TEM, 

ZIGZAG, 

PRES, 

COSMIC, 

MSCNAS, 

MHOST, 

EZRD, 

ECHO, 

RANG, DUMP, 

COMMON /DIAGS/ LDIAG(IOO) 

CONE, 

MPLOT 


COMMON /VALUES/ NP, NE, 

NDES, IMAT, KSMF, IU, 

JU, XBEG, 

XEND, 

NM0DE(3), NOPLY, 

NSECT, MESH, 

NGROUP, NCOREG, 


IGRD, MSECT(40), NC00R(3,8) 

COMMON /PLYS/ THERM0(223),MDES(223),THHF,NLEP 
COMMON /PROPTY/ EP(16),THP(5),THF(6) > AS(9) ) CS(9), 
BS(9),RA(9),RB(9),DENS 
COMMON /PROPER/ PL(71,223) 

COMMON /MFCA/ C0DES(6,40), PERT(6,40), FP(19,1,40), 

MP(15,1,40), CORC(6,5,40), DF(40), 

NF( 40), NT(2,40), NPD(40) 

COMMON /TABLES/ GNP(12,510), KEI(6,910), NCORDL(910), RAD(2,510) 

COMMON /STIFF/ EX(3,510), AB(6,5I0), DEN(510) 


This subroutine is a second level executive program. It's purpose 
is to control the generation of solid (thru-the-thickness) composite 
blades of the turboprop or turbofan types. 
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Figure 5b. - CSOLID flow chart. 
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SUBROUTINE CSREAD 


CALLED FROM ~ CSOLID 
CALLS TO — EZREAD 


ARGUMENTS (COSD, MSCD, PLYORD, HOLLOW, UNSY, PROT, EZRD, 

IMAT, NDES, NMODE, NOPLY, NSECT, NCOREG, NGROUP, 
IGRD, IU, JU, MESH, MSECT, XBEG, XEND, KREAD, KDIAG, 
KBULK, PIN, NCOOR, MAXPLY, LMAX, NGRD, MPLY, NPLY, 
KGRD, THWAL, PERT, CODES, NPM, PMOD, HOLD, LDIAG) 


This subroutine reads in all data specified in the SOLID option of 
the USER'S MANUAL. The order and format are as specified in the USER'S 
MANUAL. 

Data is stored in the appropriate arrays in common blocks for use 
in the program. 
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SUBROUTINE CTRAN 


CALLED FROM — WINPUT 
CALLS TO — none 


ARGUMENTS (X,Y,N,ANG,XA,YA) 


This subroutine is used to transform points in a two-dimensional 
coordinate system by translation and rotation to a new coordinate system. 

Input is the x and y arrays of the initial coordinates, the size of 
the arrays, the rotation angle and the x and y coordinates of the new 
origin. 

The initial x and y arrays are overwritten by the new transformed 
coordinate values. 

XPOUT = (Y(I)-YA) * SIN(ANG) + (X(I)-XA) * COS(ANG) 

YPOUT = (Y(I)-YA) * COS(ANG) - (X(I)-XA) * SIN(ANG) f 

K 
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SUBROUTINE DIAGNO 


CALLED FROM ~ OPTION 
CALLS TO — 


This subroutine interprets the user option DIAGNOSTIC and sets flags 
in the subroutines according to the diagnostic numbers given below. Each 
word of the logical common LDIAG is initially set to .FALSE, then the 
words whose number appears on the option DIAG line is set to .TRUE. This 
turns on the diagnostic print option in the corresponding subroutines. 
Presently the diagnostic output is designated to FORTRAN unit 6, the usual 
write unit. This output can be directed to other units by changing the 
value of KDIAG in subroutine CINIT. 

This feature currently exists in the following subroutines but may 
be added to others as needed by following the examples in these subrou- 
ti nes. 


CSREAD 

3 

FFREAD 

18 

SPRSET 

45 

TBGRD 

47 

TBG01 

48 

WREAD 

54 

WINPUT 

55 

WSPAR 

56 

WROTAT 

58 
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SUBROUTINE EZREAD 


CALLED FROM — CSREAD 
CALLS TO — FFREAD 


ARGUMENTS ( KBULK , KREAD , KDI AG , IMAT , NDES , NOPLY , NSECT , 
NC0REG,NGR0UP, IGRD, IU,JU,MESH,XBEG,XEND, 
NXSECT , NSPAR , NXSPAR, NYSPC , 

NMODE, MSECT , NSPDES ) 


This subroutine is a user selected option which provides for a sim- 
plified input format for program parameters as specified in the USER'S 
MANUAL. 

Parameters are input by parameter name followed by an equal sign with 
the following value being that of the named parameter. 
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SUBROUTINE FFREAO 


CALLED FROM — EZREAD OPTION 

CALLS TO — SMARTY 

ARGUMENTS (CH, LG, WHAT, LI , LABEL, I PAR, KBULK, KOUT, KDIAG, 
RE,IT,LDIAG) 


This subroutine is used to support the free-field input format al- 
lowed when reading OPTION cards or parameter values. 

When called from OPTION this subroutine, acting on one line at a 
time, first removes all leading blanks and then identifies the next four 
characters as one word and returns this word to the OPTION subroutine for 
interpretation. 

When called from EZREAD it performs the same function as for the 
OPTION subroutine and in addition it submits the character string fol- 
lowing the equal sign to subroutine SMARTY for determination of integer 
and real values in the string. These integer and real values are returned 
to EZREAD in the IT and RE arrays. 

FORTRAN unit KBULK is used as a temporary write/read file when a four 
character string is converted to one word. 
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SUBROUTINE FLIPIT 


CALLED FROM — PLYSHR SIGMAZ TBGRD 
CALLS TO — none 

ARGUMENTS (LU,XX,YY) 


This subroutine is used to rearrange the input arrays to subroutine 
INTRPL. If the arrays are not in accending order this subroutine reverses 
the order prior to submission to INTRPL and is used to reverse the order 
after INTRPL. 
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SUBROUTINE GNPBLK 


CALLED FROM — CSHELL CSOLID 
CALLS TO — none 

ARGUMENTS (WINDML,COSD,MSCD,TEM, 

NYPTS , NPAST , KBULK , NP , GNP , RAD) 


This subroutine generates NASTRAN GRID cards and TEMP cards and 
stores them in FORTRAN unit KBULK. 

The GRID and TEMP cards are in a format common to both COSMIC NASTRAN 
and MSC/NASTRAN. 
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SUBROUTINE INTRPL 


CALLED FROM — MIDWAL 
TBG01 

PLYSHR 

WIND 

SIGMAZ 

WINPUT 

TBGRD 


CALLS TO — none 





ARGUMENTS (L.X.Y.N.U.V.VDOT) 




COMMON /CUNITS/ KDES1 , 
KTTAB, 
KPR0P2, 
KPLOT 

KKEI , 

KDES2, 

KPOST, 

KBULK, 
KPR0P1 , 
KMHFOR, 

KREAD, 

KGNP, 

KREFL, 

KOUT, 

KDBANK, 

KDIAG, 


This subroutine is used throughout COBSTRAN for all calculations 
requiring the interpolation of a single-valued function and/or determi- 
nation of the first derivative of the function at each interpolated point 
(ref. 4). 

The method is based on a piecewise function for which a portion of 
the curve between a pair of given points is represented by a third-degree 
polynomial for a single-valued function. In this method the slope of the 
curve is determined at each given data point locally by the coordinates 
of five data points with the data point in question as a center point and 
two data points on each side. Each piece of the function representing a 
portion of the curve between a pair of given data points is determined 
by the coordinates of and the slopes at the points. 

When interpolation is made near the end points of the curve, two more 
points estimated at each end point are used to determine the slope of the 
curve. Thus the estimation is based on three data points, the end point 
in question and the two estimated data points. 

Arrays containing the x and y values and first derivative of each 
requested point are returned in the argument of this subroutine. 
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SUBROUTINE INVA 


CALLED FROM -- COMSA PLY PROP 

CALLS TO — none 


ARGUMENTS (X,MMM,K0UT) 


This subroutine is used to invert a matrix [A] using Gauss elimi- 
nation with maximum column pivoting. The subroutine returns the inverse 
in matrix [X] or an error message if the matrix is singular. 
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SUBROUTINE LOGO 


CALLED FROM ~ MAINPG 
CALLS TO — none 

ARGUMENTS (KOUT) 


This subroutine is used to generate the COBSTRAN logo, store comments 
on major program changes and designates credits and version number. 
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SUBROUTINE MATANG 


CALLED FROM — MSCBLK NASBLK 
CALLS TO — none 


ARGUMENTS (GA,GB,THCONN) 


This subroutine calculates the angle between the element coordinate 
system and the material coordinate system for 2-D elements used in COSMIC 
NASTRAN and MSC/NASTRAN. 

The material coordinate system is defined by the x-x or spanwise 
direction in the COBSTRAN structural coordinate system. The element co- 
ordinate system x-x axis is defined by the first two grid points of the 
element connectivity generated in subroutine TBGRD. The material angle 
is the angle between the structural x-x axis and the projection of the 
element x-x axis on the structural x-y plane. 


A1 = x value of point A 
A2 = y value of point A 
BI = x value of point B 
B2 = y value of point B 


9 = ATAN (IHf) 
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SUBROUTINE MATL 


CALLED FROM — none 
CALLS TO — none 


ARGUMENTS (WINDML, SPARS, COSD.MSCD.RINDV, HOLD, KBULK, 

NPAST , NE , I MAT , MESH , EX , AB , DEN , GNP , KE I , NUML , RAD ) 


NOT USED IN THE CURRENT VERSION OF COBSTRAN 
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SUBROUTINE MESHG 


CALLED FROM — CSOLID 
CALLS TO — none 


ARGUMENTS ( IU , JU , MESH , IGRD , NE , MSECT , KEI ) 


The blade mesh pattern and element connectivities are generated in 
this subroutine. Element numbers and related node numbers are stored in 
a common variable KEI ( I , J ) where J is the element number and I is the node 
number. 


Four user-optional mesh patterns are available; (1) triangular ele- 
ments with diagonal sides symmetric about the blade center line, 

(2) triangular elements with diagonal sides alternating chordwise, 

(3) triangular elements with three corner and three mid-side nodes, and 

(4) quadrilateral elements. 


Storage locations for these patterns are; 


(1) KEI ( I , J ) J 
I 


element number 
1,3,5 node numbers 


(2) KEI( I , J) J 
I 


element number 
1,3,5 node numbers 


( 3 ) 


KEI ( I , J ) 


J = element number 
I = 1,3,5 vertex numbers 
I = 2,4,6 mid-side numbers 


(4) KEI ( I , J ) J 
I 


element number 
1,2, 3, 4 node numbers 
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SUBROUTINE MHCONN 


CALLED FROM — MHPRE 
CALLS TO -- none 


ARGUMENTS (COSD, MSCD, TEM, PRES, NP, NE, NYTPS, NPAST, 
KBULK, MESH, GNP, KEI , NCORDL) 


This subroutine writes the element connectivities generated in pre- 
vious subroutines to the FORTRAN unit KBULK in the format required for 
the MHOST module. 

Format of output is 

I, Kl, K2, K3, K4 

where I is the element number and KI thru K4 are node numbers. 
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SUBROUTINE MHCOOR 


CALLED FROM — MHPRE 
CALLS TO — none 

ARGUMENTS (COSD, MSCD, NP, NYPTS, 
NPAST, KBULK, GNP, RAD) 


This subroutine writes the nodal coordinates generated in previous 
subroutines to the FORTRAN unit KBULK in the format rquired for the MHOST 
module. 

The first card generated is l *COOR t . 

Format of output is 

I, X, Y, Z, t 

where I is the element number, X, Y, Z are the coordinates and t is 
the nodal thickness. 
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SUBROUTINE MHPREP 


CALLED FROM — CSHELL CSOLID 

CALLS TO — MHCONN MHCOOR MHPROP 


ARGUMENTS (WINDML, SPARS, COSD, MSCD, TEM, PRES, 

NP, NE, NYPTS, NPAST, KBULK, KREAD, KPOST, 
KPR0P1, KPROP2, MESH, IMAT, HOLD, GNP, KEI, 
RAD, NCORDL) 


This subroutine reads in the users input to the MHOST module and 
rewrites this input file to FORTRAN unit KBULK including the generated 
nodal coordinates, connectivities, and nodal laminate material proper- 
ties. 


This subroutine interprets the MHOST module input file and appro- 
priately places the nodal coordinates, connectivities and material data 

Upon completion of deck generation this subroutine rewrites the en 
tire file to FORTRAN unit KREAD for use by the MHOST module. 
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SUBROUTINE MHPROP 


CALLED FROM — MHPRE 
CALLS TO — none 

ARGUMENTS (WINDML, SPARS, NP, NYPTS, 

NPAST, KBULK, KPR0P1, KPR0P2, GNP) 


This subroutine writes the nodal laminate properties generated in 
previous subroutines to FORTRAN unit KBULK in the format required for the 
MHOST module. 

Five lines of data are generated for each node; 

1) Node range I to I 

2) Membrane stiffness AS(1 thru 7), AS(9) 

3) Coupling stiffness CS(1 thru 7), CS(9) 

4) Bending stiffness BS(1 thru 7), BS(9) 

5) Thermal & density AX, AY, AXY, DENS 

where AX, AY, AXY are thermal coefficients of expansion and DENS is 
the mass density. 
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SUBROUTINE MHSHCN 


CALLED FROM — MHSHPR 
CALLS TO — none 

ARGUMENTS (SPARS, TEM, PRES, KKEI , KBULK, NCORDL, 
NXSECT, NETOT, NPTOT) 


This subroutine generates the MHOST input file for a SHELL/SPAR 
model. Elememt connectivities, distributed loads (pressure) and temper- 
ature are written to the MHOST input file. 
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SUBROUTINE MHSHPR 


CALLED FROM — CSHELL 

CALLS TO — MHSHCN SHELCR SHELPR SPARCR SPARPR 


ARGUMENTS (WINDML, SPARS , TEM , PRES , KBULK , KREAD , KPOST , KTTAB , 
KPROP1, KPR0P2, KKEI , IMAT, HOLD, KAD, NCORDL, 
NXSECT, NETOT, NPTOT, NP, NPTOT1) 


This subroutine is the main driver used to generate the MHOST input 
file. Calls to each subroutine listed creates the MHOST input file in 
FORTRAN unit KBULK. 
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SUBROUTINE MHSTR 


CALLED FROM — CPOST 


CALLS TO — none 


ARGUMENTS (WINDML, SPARS, PRTOUT, MODES, NTIMES, III, 
NP, NYPTS, NPAST, KMHFOR, KOUT, KREAD, 
KDES2, NSB, MSB) 


Read the nodal forces and moments from MHOST output file FORTRAN unit 
KMHFOR. Write these forces and moments to FORTRAN unit KDES2 and to KOUT 
if a printout is selected. 
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SUBROUTINE MIDWAL 


CALLED FROM — WREAD 
CALLS TO — INTRPL 

ARGUMENTS (PIN, L, K, I) 


For aircraft wing type structures, the input coordinates are usually 
provided on the external profile of the blade. 

This subroutine modifies the input coordinates such that the nodal 
points are moved to the mid-thickness line of the shell wall. Input 
surface coordinates are translated by half of the wall thickness along a 
normal to the surface. The slope of the surface at each point is deter- 
mined using the subroutine INTRPL. 

This subroutine may be by-passed by using the PROFILE option. 
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SUBROUTINE MSCBLK 


CALLED FROM — CSHELL CSOLID 

CALLS TO ~ MATANG RAAANG 

ARGUMENTS (WINDML, SPARS, C0SD,MSCD, PRES, TEM, RANG, 
NPAST , KBULK , NE , MESH , GNP , KEI , NCORDL) 


This subroutine generates MSC/NASTRAN bulk data cards CTRIA3, CQUAD4 
and PLOAD. Continuation cards for CTRIA3 and CQUAD4 are included with 
nodal thicknesses. 
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SUBROUTINE MSCMIN 


CALLED FROM — CSOLID 
CALLS TO — RINTRD 

ARGUMENTS ( KREAD , KOUT , I U , JU , NSECT , IGRD , NP , NE , 
NGNP,NELM,GNP, KEI ,NCORDL,RAD) 


This subroutine interprets a MSC/NASTRAN bulk data deck containing 
GRID, CTRIA3 and +CT3 cards and substitutes the corresponding values into 
arrays GNP and KEI. 

Nodal thicknesses are those represented on the connectivity contin- 
uation cards. 
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SUBROUTINE NASBLK 

CALLED FROM — CSHELL CSOLID 

CALLS TO — MATANG RAAANG 

ARGUMENTS (WINDML, SPARS , COSD , MSCD , PRES JEM , RANG , 

NPAST,KBULK,NE,MESH JMAT,GNP, KEI,NCORDL) 

This subroutine generates COSMIC/NASTRAN bulk data cards CTRIA1, 
CQUAD1 , CTRIA2, CQUAD2 and PLOAD. 

If imaT = 1 thru 4 TRIA2 or QUAD2 will be generated. 

If imaT = 5 or 6 TRIA1 or QUAD1 will be generated. 
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SUBROUTINE NASMAT 


CALLED FROM ~ CSHELL CSOLID 


CALLS TO — none 


ARGUMENTS (WINDML, SPARS, COSD,M$CD,RINDV, HOLD, 

KBULK , KPR0P1 , KPR0P2 , NPAST , NP , NE , IMAT , 
MESH , EX , DEN , GNP , KEI , NCORDL) 


Element material properties are generated in this subroutine in the 
format of the NASTRAN structural analysis program. 

Arrays AS, BS, EP, CS, RA, RB, and single value DENS are read into 
this subroutine from temporary disk storage for each node and stored in 
array SK 

The type and number of NASTRAN material cards generated for each 
element is determined by the user selected value for IMAT. 

Eight user-optional property selections are available; 

(1) anisotropic material properties based on reduced axial stiffness. 
One MAT2 card is created per element. 

(2) anisotropic material properties based on reduced bending stiffness. 
One MAT2 card is created per element. 

(3) anisotropic material properties supplied by user (same for all ele- 
ments). One MAT2 card is created per element. 

(4) isotropic material properties supplied by user (same for all ele- 
ments). One MAT1 card is created per element. 

(5) anisotropic material properties are based reduced axial stiffness 
and reduced bending stiffness. Two MAT 2 cards are created per element. 

(6) anisotropic material properties are based reduced axial stiffness, 
reduced bending stiffness and transverse shear stiffness. Three MAT 2 
cards are created per element. 
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(7) anisotropic material properties are based axial stiffness, bending 
stiffness and coupling stiffness. Three MAT 2 cards are created per ele- 
ment. 

(8) anisotropic material properties are based axial stiffness, bending 
stiffness and coupling stiffness and transverse shear stiffness. Four 
MAT2 cards are created per element. 

The [G] stress-strain matrix required for the NASTRAN MAT2 material 
card is generated for each element according to the equation; 



/NCON 

CG3k = ( & 

ky * ^k 

Membrane 

C(1) = 

1 / T 

Bendi ng 

C(2) = 

(12 * NCON 

Transverse shear 

C( 3) = 

1 / NCON 

Coupl i ng 

C( 4) = 

NCON / T 2 


where 


[A] = Nodal stiffness matrix 

i = Node at which properties exist 


NCON 


t = Thickness at node i 

NCON = Number of nodes per element 

k = Properties related to 1) membrane, 
2) bending, 3) transverse shear, 

4) coupling 
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SUBROUTINE NASMIN 


CALLED FROM — CSOLID 
CALLS TO — SMARTY 

ARGUMENTS (KREAD, KOUT, IU, JU,NSECT, IGRD,NP,NE, 
NGNP,GNP,KEI,NCORDL, RAD) 


This subroutine interprets a COSMIC/NASTRAN bulk data deck contain- 
ing GRID, CTRIA2 and PTRIA2 cards and substitutes the corresponding values 
into arrays GNP and KEI. 

Nodal thicknesses are calculated by averaging the thicknesses of all 
elements connected at a node. No weighting factors are used. 
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SUBROUTINE NASTDK 


CALLED FROM 

— CSOLID 

CSHELL 




CALLS TO — 

GNPBLK 

MSCBLK 

NASBLK 

NASMAT 

PQSBLK 


ARGUMENTS (WINDML, SPARS, COSD.MSCD, PRES, TEM, RANG, RINDV, 
MSCNAS , COSMIC , KBULK, KPR0P1 , KPR0P2 , KOUT, 

I , NXSECT , NSPAR, NPAST , NYPTS , NP , NE , IMAT , MESH , 
HOLD , EX , DEN , GNP , KEI , NCORDL , RAD ) 


This subroutine is a driver for the generation of both COSMIC and 
MSC/NASTRAN bulk data decks. The five subroutines called each generate 
specific NASTRAN bulk data card formats which are defined under each 
subrouti ne. 
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SUBROUTINE NOLAPR 


CALLED FROM — CSOLID CSHELL 

CALLS TO — OUTPUT PLY PROP 


ARGUMENTS (WINDML, PROT , PREP , MODES , COSD , MSCD , TEM .SPARS , 

SOLID. BIDE , KDES1 , KDES2, KBULK, KOUT , KPROP1 , KPROP2 , 
NGRID.NP, IMAT.THHF.NLEP.NPAST.NXSECT.NOPLY, 

NSPAR , NXSPAR, NYPTS , NPTOT, NETOT , CODES , PERT , RAD , 
PL.EX, AB.DEN, FP.MP.CORC.DF.NT.THERMO.MDES) 


This subroutine is a driver for the generation of equivalent homo- 
geneous anisotropic properties for each node from the ply layers and 
stack.-up order previously determined for each node. Each node is evalu- 
ated separately calling subroutine PLY and subroutine PROP. Nodal prop- 
erties are written to FORTRAN units KPR0P1 and KPR0P2. 
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SUBROUTINE OPTION 


CALLED FROM — MAINPG 
CALLS TO — FFREAD 


ARGUMENTS ( KBU LK , KOUT , KREAD , KDI AG 

COMMON /BOOL/ BIDE, 

DYNA, 

MINICH, 

UNSY, 

POST, 

MODES, 

PLYORD, 

RINDV, 

COSMIC, 

MSCNAS, 

RANG, 

DUMP, 


i\nnu | uuru , 

COMMON /PRINT1/ ID(80,5), IOPT(80 
COMMON /PRINT2/ NPRT(2,10) 


LDIAG) 


PRTOUT, 

WINDML, 

SOLID, 

PROFIL, 

PROT, 

PREP, 

HOLLOW, 

COSD, 

MSCD, 

TEM, 

ZIGZAG, 

PRES, 

MHOST, 

EZRD, 

ECHO, 

CONE, 

MPLOT 



30) 


This subroutine reads all user option cards and sets all logical 
controls in accordance with user choices. The “COBSTRAN OPTION SUMMARY" 
is printed from this subroutine. 

This subroutine is designed to promote a user-friendly approach to 
COBSTRAN input and to allow for ease of modification as new program 
options become available. 

All option cards with the exception of ENDOPTION are stored as log- 
ical values in COMMON block /BOOL/. Initially all values in /BOOL/ are 
set .FALSE. 
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SUBROUTINE 0UTP10 


CALLED FROM — CSOLID 


CALLS TO — none 


ARGUMENTS (MAXPLY, LMAX,NXSECT,MPLY,NPLY, KOUT) 


COMMON /BOOL/ BIDE, 

DYNA, 

PRTOUT, 

WINDML, 

SOLID, 

MINICH, 

UNSY, 

PROFIL, 

PROT, 

PREP, 

POST, 

MODES, 

HOLLOW, 

COSD, 

MSCD, 

PLYORD, 

RINDV, 

TEM, 

ZIGZAG, 

PRES, 

COSMIC, 

MSCNAS, 

MHOST, 

EZRD, 

ECHO, 

RANG, 

DUMP, 

CONE, 

MPLOT 


COMMON /PRINT1/ ID(80,5), 

IOPT(80,30) 




COMMON /PRINT2/ NPRT(2,10) 

COMMON /VALUES/ NP, NE, NDES, IMAT, KSMF, IU, JU, XBEG, XEND, 

1 NMODE(3), NOPLY, NSECT, MESH, NGROUP, NCOREG, 

2 IGRD, MSECT(40) , NC00R(3,8) 


This subroutine generates the output page titled 
"COBSTRAN INPUT SUMMARY" indicating all program options, parameter values 
and ply stacking sequence to be used by COBSTRAN. 
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SUBROUTINE 0UTP11 


CALLED FROM — CSOLID 
CALLS TO — none 

ARGUMENTS ( NDES , KOUT , CODES , FP , MP , CORC ) 


This subroutine generates the user requested print-out of the con- 
stituent fiber and matrix properties used in the analysis. This enables 
the user to examine the values stored in the COBSTRAN internal databank. 
If modifications are required to these values to meet user requirements, 
these values may be copied to an external databank, modified and used in 
the analysis. The same names may be utilized in the external databank 
as in the internal databank because COBSTRAN will select the external 
databank prior to the internal databank. This subroutine is activated 
by the DUMP option. 
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SUBROUTINE OUTPUT 


CALLED FROM — CPOST CSHELL CSOLID NOLAPR WREAD 


CALLS TO — none 


ARGUMENTS ( I I I , WI NDML , PROT , PREP , MODES , COSD , MSCD , TEM , SPARS , 
NPAST , NXSECT , NSPAR , NXSPAR , NYPTS , NPTOT , NETOT , 

KDES1 , KBULK, KOUT, KPROP1 , KPROP2 , CODES , PERT) 

COMMON /VALUES/ NP, NE, NDES, IMAT, KSMF, IU, JU, XBEG, XEND, 
NMODE(3) , NOPLY, NSECT, MESH, NGROUP, NCOREG, 
IGRD, MSECT(40) , NCOOR(3,8) 

COMMON /TABLES/ GNP(12,510),KEI(6,910),NC0RDL(910),RAD(2,510) 
COMMON /PLYS/ THERMO(223),MDES(223),THHF,NLEP 
COMMON /PROPER/ PL(71,223) 

COMMON /PRINT1/ ID(80,5), IOPT(80,30) 

COMMON /PRINT2/ NPRT(2,10) 

COMMON /PROPTY/ EP( 16 ) , THP( 5 ) , TH F( 6 ) , AS( 9 ) , CS( 9 ) , 
BS(9),RA(9),RB(9),DENS 


Various output formats are available in this subroutine for use at 
different times throughout the COBSTRAN run. Any one of nine format 
options is selected in the argument of the call statement. 


Nine different formats are printed by this subroutine. 

1. Number of nodes/elements and element types 

2. Composite types, ply geometry and distribution 

3. Table of ply stack-up order for each node 

4. Nodal elastic properties 

5. Extended ply property output for one node 

6. NASTRAN GRID cards to FORTRAN unit KBULK 

MINICH option only (Not available in Version 1.2) 

7. Blade geometry 

8. Ply failure criteria 

9. Nodal radii of curvature 
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SUBROUTINE PLTBLK 

CALLED FROM — CSOLID CSHELL 
CALLS TO — none 


ARGUMENTS (WINDML, SPARS, C0SD,MSCD, KPLOT, NYPTS, NPAST, 
NP, NE, MESH, GNP.KEI, RAD, NCORDL) 


This subroutine generates NASTRAN GRID cards and PLOTEL cards for 
the COBSTRAN model being created. This deck, written to FORTRAN unit 
KPLOT, is used for generation of an undeformed model plot. The file is 
created by using the MPLOT option. 
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SUBROUTINE PLY 


CALLED FROM — CPOST NOLAPR 

CALLS TO — INVA 

ARGUMENTS (BIDE, KOUT , THHF , NLEP , CODES , FP , MP , CORC , DF , NT , 
THERMO, MDES, PL) 


The thermoelastic properties of the single ply are generated in this 
subroutine. In addition the actual fiber and matrix volume content, the 
ply thickness, density, and the interfiber spacing are generated. The 
equations programmed to generate ply properties using constituent prop- 
erties from the resident databank are; 


kf = (1.0 - k v ) kf 


k m = (1. 0 — k v )(1.0- k f ) 
Pt = Pf ^f + ^m 



where k y and kf are read globally and denote void and fiber volume ratios, 
respectively. The subscripts f, m and £ denote fiber, matrix and ply 
property, respectively. 

The equations programmed to generate the extensional moduli and the 
thermal coefficients of expansion are; 

[E^] = CC f ,] C Ef T[C f^3 kf + [C^l^CEp^DCC^]^ 


and 


{<y} = [C f/ ] T {a f }kf + [C m ] T {« m }k m 


2-65 ( 8 - 31 - 88 ) 



COBSTRAN Theoretical/ Programmers Manual 


The arrays in the last two equations are given by; 


and 


[E 


/, f ,m 


] = 


1 

V21 

V31 

E m 

E /22 

E ^33 

Y12 

1 

Y32 

E m 

e ^22 

e ^33 

Y13 

Y23 

1 

E ^ll 

E ezz 

E <f 33 


= c “l> “ 2 ’ “3-V, f ,m 


The arrays [E^3,[Ef], and[E m ] are generated locally in the arrays ECL, 
ECF, and ECM respectively. 

The arrays [Cf/J and CC m /] and the constants (defined subsequently) 
in them are given by; 


CCf,] = 


1 

if v f 2 1 v m21 \ 

if v f 3 1 v m31 \ 

A E mll 

A V c f E f22 c m E m22 ) 

A V C f E f33 c m E m33 ) 

0 

1/C f 

0 

0 

0 

1/C f 


and 


1 

1 / v m21 v f 21 \ 

1 / 

v m31 v f 3 1 \ 

BE fll k f 

B V c m E m2 2 c f E f22 ) 

B [ 

- c m E m33 c f E f33 ) 

0 

i/c, 


0 

0 

0 


i/c. 
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where 


A = 


E fll E mll 


B = 


_JL_ + _l_k 

E mll E fll m j 



/?f =1.0 


_ a 


(7 1.0 \l;C0RC(2,l) 




"m 


) 


C0RC(2, 1) 


If C0RC(2, 1) # 0 
if C0RC(2, 1) = 0 


where the array C0RC(m,n) contains the respective correlation coeffi- 
cients available in the resident databank. 
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The variables CORC are the empirical (adjustment) factors and are 
part of the internal databank. The elements in the array CORC constitute 
experiment-theory correlation (semi-emperical ) factors and are selected 
so that the predicted and experimental results for a particular fiber- 
matrix system from a particular fabrication process are in good agreement. 
The variable Bf could be selected to be different from unity if addi- 
tional adjustment is needed. 


The elements in the arrays CC^] and [C^j are generated by substi- 
tuti ng 


Pf = 1-0 


Pm — 


1.0 U/C0RC(4.1) 
k m / 

C0RC(4, 1) 


if C0RC(4, 1) * 0 
if C0RC(4, 1) = 0 


The equations programmed for the ply heat capacity and the ply heat 
conductivities are 


H L 

M cn - p r ,• 


(H c f pf k-p + H cm p m k m ) 




^ ^kv ^maa + k v k y) 

^^m«« + k v “ k v (K maa - K v ) 


and 


K ^l = ^kl k f K fll + k m K mll 


The subscript a takes the values (1,2 and 3). The remaining variables 
are read in globally. 
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This subroutine is also used to calculate the ply conductivities 
K ^22 and K ^.33 . The specific equation programmed in this subroutine is 


1 ^ko^/ kf + 

1 _ 

^kaVk^: 

where a takes the values 2 and 3. 

The variables are fiber conductivity K faa , matrix conductivity 
K m (modified for void effects), actual fiber volume ratio kf , corre- 
lation factor /? ka , and the computed conductivity K^ aa . 


BTA = (/? kv , /? kl , /? k2 , /? k3 ) 
and 


K f 11 K f22 K f33 H cf 

CHK= K mll K m22 K m33 H cm 

0 0 0 * v 

Inputs to subroutine PLY are fiber and matrix material properties 
and the correlation factors. These properties are read in globally and 
are (E, v, Pt H c , K, a) f>m , N f , d f , VCF, BTA, and (k y , k f ) i where i=l to 
and is the number* of layers). For corresponding code identifiers 
see APPENDIX A. 

The outputs of subroutine PLY are the basic ply properties 
(kf, k m , t£ and 6 e ) which are stored in PL(3,I) and PL(5,I) to PL( 8 , 1 ) ; 
the ply stress-strain relations, which are stored in PL( 15,1) to PL( 23 , I ); 
the ply thermal coefficients of expansion, heat conductivities, and heat 
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capacity, which are stored in PL(24,I) to PL(30,I), and the ply moduli 
and Poisson’s ratios, which are stored in PL( 3 1 , I ) to PL(42,I). 

The properties for PL(1-42,I), PL(49, I) and PL(50,I) are generated 
in, or input to, this subroutine and stored in the PL array. 
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SUBROUTINE PLYIN 


CALLED FROM — COMSA CSHELL CSOLID 
CALLS TO — none 

ARGUMENTS (J,POST,KOUT,NPM,PMOD, 

CODES , FP , MP , CORC , PL) 

COMMON /PLIN/ FPT ( 19 ) , MPT (15),CC(6,5), TS3( 3 ) ,ALPP(3) 


This subroutine provides for user-intervention in the generation of 
layer p if the properties of one or more designated layers are known or 
if the effects o layer property changes are to be studied. 

Properties required for input are; fiber volume ratio,, density, 
longitudinal modulus, transverse modulus, shear modulus, poisson's ratio, 
longitudinal tensile strength, longitudinal compression strength, trans 
verse tensile strength, transverse compression strength, intralaminar 
shear longitudinal thermal expansion coefficient, and transverse thermal 
expansion coefficient. 

This subroutine modifies the properties given in the databank for 
T300 fiber and SPOX matrix so that the combination of the modified prop- 
erties when processed through the subroutines PLY and PROP will generate 
a layer with properties equivalent to the input layer properties. 
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SUBROUTINE PLYSHR 


CALLED FROM — CPOST 

CALLS TO -- FLIPIT INTRPL SIGMAZ 


ARGUMENTS ( NP , IU , JU , IGRD , PRTOUT, KDES1, KTTAB , KOUT , 

CODES , MSECT , GNP , FC) 

COMMON /PLYS/ THERMO(223),MDES(223),THHF,NLEP 
COMMON /STRA/ DISV1(6,510) 

COMMON /PRINT1/ ID(80,5), IOPT(80,30) 

COMMON /PRINT2/ NPRT(2,10) 

COMMON /STIFF/ VF(2,510), XINT(36), YINT(36), 

XOUT(36) , XM( 14, 36) ,DUM( 14) , VFP(2, 510) , 
SHL(8, 223) , SCXZ(223) , SCYZ( 223) , DUMM1(204) 


Interply shear and tensile stresses are calculated in this subrou- 
tine. Also interply limiting shear strengths and margins of safety are 
computed. 

The nodal shear forces V x and V y are determined by differentiating 
nodal bending moments using subroutine INTRPL. 

The sign convention used for this analysis is shown in Figure. 6a. 


Summing moments yields; 



3M X y 

ay 


Qy 


aM y aM xy 
"ay ax - 
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The interply shear stresses are calculated according to the 
parabolic shear stress distribution for rectangular cross sections. 


-l[(h) 2 -y§] 

where I = b( | ^ - b = 1 

12 

therefore 



The relationship of b, h and yp is shown in figure 6b. 
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Figure 6b. - Rectangular cross-section variables. 
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The interply tensile stresses are determined as follows; 
Force equilibrium equations 


da x ^y_ d(T *y = o 

d X dy da 2 


da 


yx 


dx 


+ 


d<Ty 

ay 


da 


+ 


yz 


dz 


= 0 


da zx 

dX 


+ 


da zy 

dy 



= 0 


rewriting the last equation 

da z _ da zx d<7 zy 

dz dx dy 


where 



and where the shears Q x and Q y were defined previously, 
therefore 



aQx. + J^y 

dX dy 


3_ J_ 

4 h 


1 - 



integration of the above equation yields 
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d z = 



yo.__i 

h 3 



+ C 


Delamination Limiting Strength 

When the composite is subjected to loads normal to its plane, or have 
rapid moment changes, interply layer delamination could be caused by the 
transverse shear stresses ^cxz anc * a cyz* 

The limiting condition on these stresses is established by the pro- 
cedure discussed in Ref. 2, Chapter 2 for obtaining the limiting condition 
on the ply's intralaminar shear stresses. Adopting that procedure, it 
can be shown that 


where 


S ._ B f mpS 

5 cxzj-Pcxzl ^ cxz rexz 




E mps 


'cyzsj-^cyz\ 4>^ cyz ^ cyzj 


0 


^CXZ 


Vi + 1 + Vi 

r r (Vi + 1 “ Vi+l^mi+l , (Vi “ Vi> G mi 

Vi+1 + Vi + ~ + 


Vxzi + 1 


Yxzi 


Vcyz 


Vi+1 + Vi 


f , e . (Vi + 1 Vi+l^mi + 1 , (Vi“Vi^mi 

Vi+1 + Vi + ~ + 


Yyzi+1 


Yyzi 
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The moduli G cxz j and G C y Z j are the local transverse shear moduli in 
their respective planes, from (Ref. 2). 


G 


cxzj 


Vi+1 + Vi 
Vi + 1 Vi 

Vxzi+1 Vxzi 


G 

u cyzj ~ 


Vi + 1 + Vi 
Vi+1 Vi 

Vyzi+1 Vyzi 


The parameters P cxz and /? cyz are the theory-experiment correlation 
factor. Its numerical value is determined from horizontal shear beam 
tests on composites with coincident material and structural axes 
(Ref. 2). 

The shear moduli G^ xzl - and G^.y Z1 * are the ply transverse shear moduli 
G^j 3 andG ^23 transformed along the structural axes. 
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SUBROUTINE PQSBLK 


CALLED FROM — NASTDK 
CALLS TO ~ none 


ARGUMENTS (WINDML, SPARS, COSMIC, C0SD,MSCD, 

KBULK, KOUT, NPAST, NE, IMAT, MESH, 
GNP,KEI,NCORDL) 


This subroutine generates NASTRAN property (PSHELL cards) for 
MSC/NASTRAN QUAD4 and TRIA3 elements. Factors on the PSHELL card are left 
blank. MSC/NASTRAN default values will be used in the analysis unless 
these cards are edited by the user. 

For COSMIC/NASTRAN PQUAD1, PQUAD2 , PTRIA1 and PTRIA2 property cards 
are generated for the corresponding element cards. Factors on the P cards 
are; 


T1 = t, T2 = (t 3 )/12, T3 = ( 2/3 ) * t 


Material property ID numbers are generated and inserted on each 
P card according to the users selected parameter IMAT. 

Material property ID's are; 

Series 

100000 Membrane or Reduced Membrane 

200000 Bending or Reduced Bending 

300000 Tranverse shear 

400000 Coupling 
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SUBROUTINE PROP 


CALLED FROM — 

CPOST 

NOLAPR 

CALLS TO — 

INVA 



ARGUMENTS (II, KOUT, KPR0P1 , KPR0P2 , NP , IMAT, THHF, NLEP , 
NGRID, NYPTS , NPTOT , P L , EX , AB , DEN ) 

COMMON /PROPTY/ EP(16),THP(5),THF(6),AS(9),C$(9), 
BS(9) , RA(9) , RB(9) ,DENS 


This subroutine generates the three-dimensional thermoelastic prop- 
erties of the composite about its structural (x,y,z) and material (1,2,3) 
axes. The angle 0 is measured from the x axis of the structural axes 
system to the material axis of each ply (See Figure 7). These composite 
properties are generated from the following equations; 

N e N,-l 

/^N z .fi+1 ~ ] + 2] NjCSj] 

i=l J=1 

N, 

(“cl = “t [E c ] ^ ^ ( z ^1+l ~ 
c i = l 

The arrays {« c } and in the preceeding equations are given by 

{a c } = t« CX x a cyy a czz a cyz a czx a cxy^ T 



and 

{«^} = [ a ^H <y22 V33 0 0 0]T 
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Figure 7. - Ply orientation geometry. Composite structural axes, x, y, z; 
composite material axes, 1 , 2, 3; ply material axes (coincides with fiber 
direction), 1| , 2| , 3 ( . 
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For all practical purposes the two-dimensional thermal coefficients 
of expansion about the composite structural axes are the same as 
« cxx . a cyy» “cxy in the array for the three-dimensional case. 


The matrices CE c ], [E n O, [R n -1, and [Sj] are given by 



1 

v cZl 

v c31 

E cll 

E c22 

e c33 

"cl 2 

1 

v c32 

E cll 

e c22 

e c33 

v cl3 

v c23 

1 

E cll 

E c22 

E c33 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 


0 



0 


0 

0 

0 


0 


1 

e c31 

0 
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Note that for the case of an anisotropic material, the elements (1,6) 
(2,6) (3,6) (4,5), and their symmetric parts will not be zero. 


[E n O 


-1 


1 Y21 

E m £ az 

VI 2 1 

E m e ^22 

VI 3 V23 

E eZZ 

0 0 

0 0 

0 0 


o 

*Y33 

Y32 o 
E /33 


E /33 



0 0 
0 0 


0 0 

0 0 

0 0 

0 0 
1 

E ,31 



[R,i] 


cos ^0 sin^O 0 0 

sin^0 cos^o 0 0 

0 0 10 

0 0 0 cos o 

0 0 0 - sin 0 

- s i n Zo s i n Zo 0 0 


0 1/2 si n Zo 

0 -1/2 sin Zo 

0 0 

sinfl 0 

cos o 0 

0 cos Zo 


where 0 = 0 ■** for properties about the composite material axes and 
0 = 0^} +0 CS for properties about the composite structural axes (see 
Figure 7). 
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[SJ = 


( sin - sin - ( sin 20 j - sin 2o^_^ 000 - ( sin 20 - • - sin 2o. |_j) 

x( cos - cos 2^* 

( si n 20^ - si n 2 (). *_i) 
x( cos 2tf.j - cos 2fl •,*_]_) 


( sin 20^ - sin 2o^_^ ( sin 20^ - sin 2(? 1 *_^)^ 000 ( sin 20j - sin 2^-j ^) 


000 

000 

000 


- ( sin 2 0} - sin 20- -_j) ( sin 2fl.j ~ sin 2fl.j_i) 000 ( cos 20-j - cos 2fl-j__i) < 
x( cos 2tf.j - cos 2dj_i) x( cos 20 j - cos 20^^) 


where i>l and denotes the ply index. The angles (). • and0-j_j (Figure 7) 
are given by 


^i ~ *Vi + ^cs 
°i-l = *Vi-l + fl cs 
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The composite heat capacity is the same for both the three- and 
two-dimensional cases. It is given by 



Z tviVi 

i=l 


and t c is given by 


t c = 'Tj Vi 
i = l 

The composite three-dimensional heat conductivities along the com- 
posite material axes, assuming an orthotropic composite, are given by 


V 

K c n = J- cos2( V + hn s1 " 2( V)i 

° i=l 

Nf 

K C 22 = T" sin2( V + co$2 ^)i 

C 1=1 

i = l 

.Jhe angle 0 # measured from the material axes (fig. 7). 
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The composite two-dimensional heat conductivities along the compos- 
ite structural axes are given by; 


K cxx =T~ E ViOVll co%2 ° + *r2Z sin2 «)i 
c 1=1 

N, 

K cyy = -f- E Vl'( K /ll sin2fl + K /22 cos2fl )i 
c i = 1 

N, 

K cyx = K cxy = ~t~ S Vi ^22 " VllM sl * n 2 ^i 
c i = 1 

^ czz = ^c33 

The angle 0 in the last set of equations is measured from the com- 
posite structural axes and is equal to 0 cs + 0^ . 

This subroutine generates the required section properties and the 
force-deformation-temperature relations for a two-dimensional multilay- 
ered composite. It also generates the plane-stress elastic constants for 
the composite. The force-deformation-temperature relations generated in 
this procedure are defined in the following equation; 



[A C x-^ J_ ^ f £ csx i f jjjcATxj 1 

[C cx ] i CD cx ] J l w cbx j | {M caTx } J 
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The generic equations for the elements in the arrays 
CAcxl. * {N cA j x }, an< ^ ( m caTx> are 


N/ N,-l 

^cx-^ = X A T^i( z jfi+l ~ z ^i)[R/i^CE^>i] HjCSj] 

1=1 j=l 


N, N,-l 

^cx-* = ~"^~ X# A T^i(zfi + i - ^ z rpj^j^j-* 

i=l j=l 


CDcx] = -| 


N/ 

5 ] A ^i( 2 /i + l - ^CR^j] + 

i=l 


N,-l 

Yj Z ?pj H j CS j ] 
j=l 


{^caTx} = AT^i(z^i + j 
1=1 


-z n .)CR n .][E n .] 


1 


{<Vi> 


{^caTx} = ~ \ X A ^Vi( z /1+l ” z lii^i^[fviJ ^{<y 0 
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The arrays (a^-j), C R^-j H , CE^j] , and [Sj] are; 


7T 

II 

i — i 
R 

11 “22 0] i 



2 

COS 0 

sin 2 0 

1/2 si n 2o 


. 2 
sin^tf 

? - 
cos 0 

1/2 sin 20 

- s i n 2o 

sin 2 0 

cos 2 o 


1 

V21 

0 



E m 

E eZZ 



V12 

1 

0 



E «?1X 

e <?22 



0 

0 

1 



<V12 

i 


S j 22 = s jll = -J- ( sin 2o i - sin 2 /? i-l ) 2 
s j21 = s jl2 = - s jll 


s j32 = S j23 = x ( 5111 2o i “ s1n 2o i — 1)( cos 2w i “ cos 2fl i-l) 


s j31 = s jl3 = - s j23 



( cos 20 j - cos 


2 


Here 0-j equals the o cs + Qj in figure 5. 

The reduced bending rigidities (ref. 5) are generated in this pro 
cedure according to the equation 
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The reduced axial stiffnesses are generated in the procedure ac- 
cording to the equation 




'CX 




The two-dimensional composite elastic constants are generated from 
the following equation (assuming aT^- = aT for i=l to N^ ); 


[E 

LC cx J ~ + 


Z<Vi+l 

i=l 




IR/i 1 + 


2 

j = l 


¥ S j ] 


where 


■V 

Vi 

i = l 

The inputs to this subroutine are N^, + j , z^ , 0 CS , 0^ , CE n -], 

Hj, {a^} , h^j , and {K^.,-} , which are all global. Other inputs are t^- , 
AT n - , and dj (relative to the composite structural axes), and the ply 
elastic constants. These quantities are global and are located, respec- 
tively in PL(7, I) , PL(50, I) , PL(14,I), PL(9,I), and PL(31,I) to PL(42,I). 
The remaining quantities are either generated or are transferred from 
information stored in PL(8,I), PL(11,I), P L( 13 , 1 ) , PL( 15-23 , 1 ) , 

PL(24-30, I) . 

The arrays [R^.|]T, [E^] - * , [R^-j], and [Sj] and the dimensions z n - 
are generated in this subroutine. 


The outputs are t c and the arrays CE c ]~l, {« c }, [E C T, h c , and {K c }. 
The composite thickness t c is stored in GNP(6,I). The arrays 
[E c ]~l, and {a c }, for the composite structural axes are printed out under 
the heading; 

ELASTIC PROPERTIES — (STRAIN-STRESS TEMPERATURE RELATIONS) 
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The composite structural axes properties [E c ]~l and {« c } are stored 
in EP(1) thru EP( 16) as global variables. The heat capacity and three- 
dimensional heat conductivities along the structural axes are stored in 
THP(l) thru THP(5). The two-dimensional thermal coefficients of expan- 
sion along the structural axes are stored in EP( 13) , EP( 14) and EP(16). 
The two-dimensional heat conductivities and heat capacity along the 
structural axes are stored in THP(l) thru THP(4). Note that the heat 
capacity is a scalar quantity and is independent of the reference axes. 

The density, heat capacity and three-dimensional heat conductivities 
stored in array THP are printed under the heading; 

THERMAL PROPERTIES 

The axial stiffnesses stored in array AS are printed under the 
heading; 


AXIAL STIFFNESS 

The bending rigidities stored in array BS are printed under the 
headi ng; 


BENDING STIFFNESS 


The coupling stiffnesses stored in array CS are printed under the 
heading; 


COUPLING STIFFNESS 


The reduced axial stiffnesses stored in array RA are printed under 
the heading; 


REDUCED AXIAL STIFFNESS 

The reduced bending stiffnesses stored in array RB are printed under 
the heading; 


REDUCED BENDING STIFFNESS 


The outputs are the force-deformation-temperature relations, which 
are stored in the global arrays AS=A CX , RA=A§ X , CS = C cx , 

BS = D cx , RB-D^ X , NSDT-N cA j x . 

The distances L,- and z ^ , are stored, respectively in PL( 10 , 1 ) , and 
PL(ll.I). 
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SUBROUTINE RAAANG 


CALLED FROM — MSCBLK NASBLK 
CALLS TO — none 


ARGUMENTS (GA,GB,GC, RAANG) 


The purpose of this subroutine is to calculate the angle between the 
element coordinate system and the projection of the COBSTRAN structural 
x-axis on the plane in which the element lies for 2-D elements used in 
COSMIC NASTRAN and MSC/NASTRAN. 

The material coordinate system is defined by the x-x or spanwise 
direction in the COBSTRAN structural coordinate system. The element co- 
ordinate system x-x axis is defined by the first two grid points of the 
element connectivity generated in subroutine TBGRD. The material angle 
is the angle between the projection of the structural x-x axis onto the 
plane of the element and the element x-x axis (fig. 8). 


A -> B element x-axis 

A -> N projection of the structural x-axis on the 
plane in which the element lies 

Solve for the angle 0 

o is positive if N is on the C side 
0 is negative if N is opposite C side 

Let the coordinate of each point be represented by the tag 1,2,3 for 
x,y,z respectively. 


therefore; 

A1 = x value of point A from the local origin 
B2 = y value of point B from the local origin 
etc. 

The origin lies in a plane normal to the structural x-axis and con- 
tains points B and C. 

Solving for the coordinates of N; 

N1 = B1 = Cl = 0.0 

Slope of C-B line is m = ■■■ £§— - I I- 

C3 — B3 
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A1 so 


N2 - 

B2 = m(N3 -B3) 

(i) 

N2 = 

-4-N3 

ID 


(2) 

N3 = 

- m(N2) 


(3) 

Substituting equation 

(3) into equation (1) 


N2 - 

B2 = m(-mN2 • 

-B2) 


N2 = 

-mB3 + B2 




1 + m2 



N3 = 

- mN2 




Let AA, BB, and NN represent the length of the lines opposite the 
A, B, N nodes. 


AA = V(B1 - Nl)2 + (B2 - N2)2 + ( B3 - N3)2 


BB = v(Al - Nl)2 + (A2 - N2)2 + (A3 - N3)2 


NN = V(A1 - Bl)2 + (A2 - B2)2 + (A3 - B3)2 


SS 


AA + BB -f NN 

2.0 


R = 



° = 2 -°( atan ^aa) 
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X 



Figure 8. - Determination of material x-axis relative to element x-axis. 
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SUBROUTINE REFLEC 


CALLED FROM — MAINPG 
CALLS TO — none 

ARGUMENTS (N, KREAD, KREFL.KOUT) 


This subroutine stores the users input data set in FORTRAN unit 
KREFL. If subsequent interpretation of the input, indicated by the ECHO 
option, calls for the printing of the user input data set this subroutine 
is again called for printing. 
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SUBROUTINE RINTRO 

CALLED FROM — MSCMIN NASMIN 
CALLS TO — none 

ARGUMENTS (L,M,N,R) 

NOT USED IN THE CURRENT VERSION OF COBSTRAN 
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SUBROUTINE SHELCR 


CALLED FROM — MHSHPR 
CALLS TO — none 

ARGUMENTS ( SPARS , NP , NYPTS , NPAST , GNP , KTTAB , KBULK) 


This subroutine writes coordinates for shell grid points to FORTRAN 
unit KBULK and coordinates for spar grid points to FORTRAN unit KTTAB. 
Grid point coordinate data is written in the MHOST input format. 
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SUBROUTINE SHELPR 

CALLED FROM — MHSHPR 
CALLS TO — none 

ARGUMENTS (NP, NYPTS, NPAST, KBULK, KPROP1, KPR0P2, GNP) 

This subroutine writes the membrane, bending and coupling properties 
of each shell grid point in the MHOST format to FORTRAN unit KBULK. 
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SUBROUTINE SIGMAZ 


CALLED FROM — PLYSHR 

CALLS TO — FLIPIT INTRPL 


ARGUMENTS ( IU,MSECT,GNP) 

COMMON /STIFF/ VF(2,510), XINT(36), YINT(36), 

X0UT(36) , XM( 14, 36) ,DUM( 14) , VFP( 2 ,510) , DUMMY(2434) 


Derivatives of the nodal shear forces in the spanwise and chordwise 
direction are calculated in this subroutine. These nodal shear forces 
are differentiated with respect to x and y using subroutine INTRPL. 

Determination of the derivatives of the shear forces is required to 
calculate the tensile stresses thru-the-thickness. 

The shear derivatives are returned to subroutine PLYSHR to calculate 
thru-the-thickness tensile stresses. 
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SUBROUTINE SMARTY 


CALLED FROM — FFREAD MSCMIN NASMIN 
CALLS TO — none 


ARGUMENTS (L,M,NF,RF,NL,KOUT) 


This subroutine is used to interpret a string of integer or real 
numbers in a free-field format. The values may be separated by commas 
or blanks. Integer values are returned in array IT and real values are 
returned in array RE. 

This subroutine is used to interpret a given field of characters in 
which a real or integer number is stored. Integer numbers do not have to 
be right justified. Trailing blanks are not interpreted as zeros. 

This subroutine is used specifically to evaluate eight character 
fields on NASTRAN bulk data cards. 
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SUBROUTINE SOLVEQ 


CALLED FROM — MAINPG 

CALLS TO — HOST, BLOCK, INTINT 


ARGUMENTS ( ISIZE, IWORK) 

C COMMON / / IWORK ( 100 000 ) 

COMMON / MACHIN / IDP 
DIMENSION IWORK ( 1 ) 

C DATA ISIZE / 100 000 / 

VERSNO = 4.2D0 


This subroutine is main driver for the MHOST general purpose finite 
element structural analysis code. 
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SUBROUTINE SPARCR 


CALLED FROM — MHSHPR 
CALLS TO — none 

ARGUMENTS (NP, NPTOT, GNP, KAD, KTTAB, KBULK) 


The local spar node number in GNP(9,I) is incremented by the total 
number of shell nodes for the creation of duplicate node numbers as re- 
quired input to MHOST. 

IREC is the shell node number corresponding to the spar node. 

Write the new spar node number, coordinates and thickness to FORTRAN 
unit KBULK for generation of duplicate nodes. 
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SUBROUTINE SPARPR 


CALLED FROM — MHSHPR 
CALLS TO — none 

ARGUMENTS (NP, NPTOT, KBULK, KPR0P1, KPR0P2, GNP) 


This subroutine writes the membrane, bending and coupling properties 
of each spar node to FORTRAN unit KBULK in the MHOST format. 


2-101 (8-31-88) 



COBSTRAN Thcoretical/Programmcrs Manual 


SUBROUTINE SPRSET 


CALLED FROM — WINPUT 
CALLS TO — none 


ARGUMENTS ( I , K, KK, NNQ, KOUT, KDIAG, NS PAR, DELT, 
ZIGZAG, SY, YINT, LSP, LDIAG) 


Thi s 
chordwi se 


subroutine determines the spar positions and modifies the 
grid positions to satisfy the user requested spar positions. 


2-102 ( 8 - 31 - 88 ) 



COBSTRAN Theoretical/ Programmers Manual 


SUBROUTINE SSTRAN 


CALLED FROM — STR 
CALLS TO — none 


ARGUMENTS ( X , Y , N , ANG , KDI AG , LDI AG) 


The purpose of this subroutine is to transform two-dimensional 
stress relations to a new coordinate system. 

The stress transformation matrix is generated in this subroutine and 
element stresses are transformed from the element coordinate system to 
the COBSTRAN structural coordinate system. 

The use of option DIAG = 49, in post processing, will print the 
transformation matrix for each element of the analysis. 



2 

cos ^0 

. 2 
sin^e 

- sin Zo 

sin 2 0 

2 

COS^O 

sin Zo 

sin Z 0 

- 4- sin Zo 

cos Zo 


2 
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SUBROUTINE $TR 


CALLED FROM — CPOST 
CALLS TO — SSTRAN 

ARGUMENTS (NTIMES, III, DYNA,PRTOUT,WINDML, MODES, RINDV,-KPOST, 
KOUT , KDES2 , NP , NE , MESH , NMODE , NPAST , NXSECT , NYPTS ) 

COMMON /STRA/ DISV1(6,510) 

COMMON /TABLES/ GNP(12, 510) , KEI(6, 910) , NC0RDL(9I0) , RAD(2, 510) 


This subroutine reads element stresses from a NASTRAN punch output 
file and calculates the applied membrane loads and applied bending moments 
at each node. 

The punch format of NASTRAN is used to store and recover the element 
two-dimensional stresses at each surface. When read by this subroutine 
the stresses on the positive z surface are stored in array XS1 and the 
stresses from the negative z surface are stored in array XS2. 

Subroutine SSTRAN is called to transform these stresses on both 
surfaces from the NASTRAN element coordinate system to the COBSTRAN 
structural coordinate system. 

The applied membrane loads and applied bending moments at each node 
are determined from the following equations; 



k 

{MSB} = ^ ^ t j 1 j - o-2-j) 

1=1 
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where 

k = Number of elements at node 
t = average element thickness 

crl, aZ are 2D stresses on upper and lower surface, respectively. 

The inputs to this subroutine are element thicknesses and two- 
dimensional stresses at each surface. 

The output is two-dimensional nodal membrane loads and bending mo- 
ments stored in arrays NSB and MSB and printed out under the heading; 


NODAL FORCES AND MOMENTS PER UNIT WIDTH IN STRUCTURAL COORDINATE SYSTEM 

^cxx> ^cyy» ^cxy-^ 

^cxx» ^cyy* ^cxy-^ 
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SUBROUTINE TBG01 


CALLED FROM — TBGRD 
CALLS TO — INTRPL 


ARGUMENTS (NP, NE, IU, JU, MESH, XBEG, XEND, NSECT, 
NPOS, INFLAG, MSECT, PIN, TIN, XINT, YINT, 
XOUT, YOUT, VDOT, KDIAG, LDIAG) 


This subroutine is utilized if IGRD = 1. For that case all chordwise 
grid points are divided into JU-1 equal spaces. All other grid point re- 
lationships such as z coordinate, thickness and temperature are interpo- 
lated. 

The blade is then divided equally spanwise into IU- 1 spaces and all 
other grid point relationships are interpolated with respect to the x-axis 
position. 

IU and JU are user input parameters. IU is the number of cross- 
sections spanwise and JU is the number of grid points chordwise. 
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SUBROUTINE TBGRD 


CALLED FROM — 

CSOLID 



CALLS TO — 

FLIPIT 

INTRPL 

TBG01 


ARGUMENTS (TEM, PRES, NP, NE, IU, JU, MESH, 

XBEG, XEND, NSECT, NCOREG, IGRD, MSECT, PIN, 
GNP, NCORDL, NCOOR, RAD, KDIAG, LDIAG) 


This subroutine reads blade profile data including temperatures and 
pressures on both surfaces. The usual input form for the blade profile 
is x (radial coordinate), y (chordwise coordinate) and, at each y value, 
two values of z, one for each surface. 

These input values may be used directly as nodal properties over 
which a mesh will be created (IGRD=2) or this input may be used as a bases 
for interpolating or extrapolating other points and loads to generate a 
finer mesh (IGRD=1). 

For IGRD=1 and IGRD=2 two values are required as input for z and the 
blade thickness is calculated by this subroutine. 

If all coordinates are to be used directly, including thicknesses 
at the nodes, IGRD=3 is to be selected. 

For IGRD = 4 the input form is the same as IGRD = 3 except that the 
feature of interpolation and extrapolation over the surface is utilized 
as in IGRD = 1. 
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SUBROUTINE TBPLY 


CALLED FROM — CSOLID 
CALLS TO — none 

ARGUMENTS ( PLYORD,UNSY, HOLLOW, COSD, MSCD, KDES1, KOUT,NDES,NP, 
IU,JU, IGRD,MSECT,MAXPLY, LMAX,NGRD, CODES, PERT, 

MP LY , NPLY, KGRD , THWAL , GNP , RAD) 

COMMON /PLYS/ THERMO(223) ,MDES(223) ,THHF,NLEP 


This subroutine determines the number and type of plys of each nodal 
point in accordance with the following input procedures; 

If the SOLID option is chosen the maximum number of plys for each 
material is based on a percentage of the thickest node or, if PLYORDER 
is specified, is based on the node thickness and plys needed to fill that 
thickness. For all nodes the maximum possible number of plys of each 
material is used starting at the surface until the node thickness is at- 
tained. As each ply is added the remaining thickness is checked until 
it is less than the thickness of the next ply to be included. If the 
remaining thickness is greater-than or equal-to the half-thickness of the 
next ply that ply is added to the analysis at that node. 

If PLYORDER option is designated, the ply order thru the second 
half-thickness will be a mirror image of the first half-thickness ply 
order specified with the first ply designated placed on the outer surface. 

If UNSYMMETRICAL option is designated, a second ply order is required 
for the opposite half of the blade starting at the outer surface. 

If the hollow option is chosen the same procedure is followed as for 
the solid option except at selected grid points the wall thickness is used 
and the remaining thickness is composed of the last ply material listed 
in the input deck. 
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SUBROUTINE WIND 


CALLED FROM — CSHELL 
CALLS TO — INTRPL 

ARGUMENTS ( I , KTTAB, NP , NE, NPAST , NNQ, NGRID, NELT , NYPTS , 
NXSPC, L$ECT,X,GNP,KEI) 

COMMON /PROPER/ TTAB(30,8,30,2),DUMMY(1433) 


The mesh pattern and element connectivities for quadrilateral ele- 
ments are defined in this subroutine. Each section is divided spanwise 
into NXSPC spaces along the x-axis and coordinate values are defined by 
interpolation between stations. 
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SUBROUTINE WINOUT 


CALLED FROM — CSHELL 


CALLS TO — none 


ARGUMENTS (SPARS , KOUT , KDES 1 , KDES2 , NP , NE , IMAT) 

COMMON /TABLES/ GNP(12,510),KEI(6,910),NCORDL(910),RAD(2,510) 
COMMON /PLYS/ THERM0(223) ,MDES(223) ,THHF,NLEP 
COMMON /PRINT1/ ID(80,5), IOPT(80,30) 

COMMON /PRINT2/ NPRT(2,10) 


This subroutine prints the element identification table, nodal ge- 
ometry table and surface data. 
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SUBROUTINE WINPLY 

CALLED FROM — CSHELL 
CALLS TO — none 

ARGUMENTS ( SPARS , PLYORD , KOUT , KDES2 , NP , NDES , THMAX , NSPAR , 
MAXPLY,NSPDES,MPLY, CODES, PERT, GNP) 

COMMON /PLYS/ THERMO(223) ,MDES(223) ,THHF,NLEP 


This subroutine determines the number and type of plys at each nodal 
point in accordance with the input procedure specified. 

The maximum number of plys for each material is based on a percentage 
of the thickest node or, if the PLYORDER is specified, is based on the 
node thickness and plys needed to fill that thickness. For other nodes 
the maximum number of plys of each material is used starting at the sur- 
face until the thickness is attained. 

The ply order thru the thickness will be a mirror image of the 
half-thickness order specified with the first ply designated placed on 
the outer surface. This subroutine is analogous to subroutine TBPLY. 
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SUBROUTINE WINPUT 


CALLED FROM — CSHELL 

CALLS TO — CTRAN INTRPL SPRSET 


ARGUMENTS ( KTTAB , KOUT , KD I AG , P I N , TTAB , S Y , MWSCT , X , LS P , 
LSECT,NGRID,NELT,NXSECT,NSPAR,NXSPAR,NYSPC, 
NYPTS , NNQ , NNX , LDI AG , ZIGZAG) 


Each input section profile is rotated and translated to a coordinate 
system with the origin at the leading edge and one axis along a chord line 
from the leading edge to the trailing edge, each section is then divided 
into NYSPC equal spaces. The space at the leading edge is then divided 
in half. 

If spars exist at the input section the spar positioning subroutine 
SPRSET is called to rearrange the equally spaced nodes so the nearest 
nodal points fall at the spar positions. Holding the spar position, the 
spaces are again equally divided between spars. 
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SUBROUTINE WREAO 


CALLED FROM — 

CSHELL 



CALLS TO — 

MIDWAL 

0UTP10 

OUTPUT 


ARGUMENTS ( WINDML , PROFI L , PROT , PREP , MODES , COSD , MSCD , 

CONE , P LYORD , TEM , EZRD , ZIGZAG , KREAD , KOUT , KOI AG , 
P I N , NPM , PMOD , CODES , P ERT , LD I AG ) 


COMMON /VALUES/ NP, NE, NDES, IMAT, KSMF, IU, JU, XBEG, XEND, 


NMOOE(3) , NOPLY, NSECT, MESH, NGROUP, NCOREG, 
IGRD, MSECT(40), NCOOR(3,8) 

COMMON /WINDS/ $Y(4,30,5), 

MWSCT(30,2) 

, ALPHA(3, 3) , 

X ( 30 ) , 

NXSPC(30) , 

NSPDES(5), 

NGRID, 

NELT, 

NPAST, 

THMAX, 

NXSECT, 

NSPAR, 

NXSPAR, 

NYSPC, 

LSP(2,5,30), 

NYPTS, 

NNQ, 

NNX, 

LSECT(2) , 

NPTOT, 

NETOT, 

HOLD(IO), 

MAX PLY, 

MPLY(300) 


This subroutine reads in all data specified in the SHELL/SPAR option 
of the USER'S MANUAL. The order and format are as specified in the USER'S 
MANUAL. 

Data is stored in the appropriate arrays in common blocks for use 
in the program. 
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SUBROUTINE WROTAT 


CALLED FROM — CSHELL 
CALLS TO — none 

ARGUMENTS ( KK , NP , KDI AG , ALPHA, GNP , LDIAG) 


If CONE is designated in the user options this subroutine calculates 
a coordinate transformation from a blade oriented coordinate system to a 
rotational axis coordinate system for all nodal position values using the 
transformation matrix ALPHA provided by the user. 
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SUBROUTINE WSPAR 


CALLED FROM — CSHELL 
CALLS TO — none 


ARGUMENTS ( KDIAG , NP , NE , NELT , NSPAR, NYPTS , NXSPAR , 
X , LSP , NXSPC , SY , GNP , KEI , LDIAG) 


This subroutine generates the connectivity and grid points for spar 
elements utilizing the existing shell grid points. Spar thicknesses are 
also evaluated in this subroutine. 
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SYMBOLS and CODE IDENTIFIERS 


Engi neeri ng 
symbol 

FORTRAN code 

Comment 

\x 

AS 

composite axial stiffness 

a r 

A cx 

RA 

reduced axial stiffness 

BIDE 

Boolean 

TRUE if interply effects 
are included 

C cx 

CS 

composite coupling 
stiffness 

Dcx 

B$ 

composite flexural rigid- 
ities 

dr 

u cx 

RB 

reduced bending rigidities 

D v 

DISV, DISVI 

displacement vectors 

d f 

DIAF 

filament equivalent diam- 
eter 

E c 

EC I 

elastic property matrix 

E f > E cf 

ECF 

filament elastic constants 

E /> E c/ 

ECL 

ply elastic constants 

E E 
^cm 

ECM 

matrix elastic constants 

E fll> E /ll» E mn 

EFlljELll, EM 1 1 

filament, ply and matrix 
normal moduli 
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Eng i neeri ng 
symbol 

FORTRAN code 

Comment 

E fl2- E /12> E ml2 

EF12,EL12,EM12 

filament, ply and matrix 
shear modul i 

E n- 

ELS 

ply elastic matrix 

H j ' 

PL(9,I) 

interply distortion energy 
coefficient 

^kc 

CHK 

array of constituent heat 
conductivities 

h c 

HHC 

composite heat capacity 
stored in PC(18) and 
PC(54) 

i J 

I,J 

index, generally ply or 
interply 

K cll> K c22’ K c33 

HK11,HK22,HK33 

composite three- 
dimensional heat 
conductivities along the 
material axes in PC(15 to 
17) 

^cxy>^cyy’ ^cyx 

HK11 ,HK22,HK33 

composite two-dimensional 
heat conductivities in 
PC(51 to 53) 

K f 11 5 ^Vll > K mll 

CHK 

see H kc 

L-f * 

KF, KV 

apparent fiber and void 
volume ratios 


KFB, KMB 

actual fiber and matrix 
volume ratios 
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igi neeri ng 
symbol 

FORTRAN code 

Comment 

^f/» 

KFL, KVL 

ply apparent fiber and 
void volume ratios 

*-sc 

LSC 

array of limiting condi- 
tions 

^caTx 

MSDT 

thermal moments 

^cx 

MSB 

applied moment 

m 

M 

load condition index 

^cx 

NSB 

applied membrane loads 

^caTx 

NSDT 

thermal force 

N f 

NFPE 

number of filaments per 
end 

N/ 

NL 

number of plies 

N/c 

NLC 

number of load conditions 

Npc 

NPC 

string PROPC length 

Npl 

NPL 

string PROP length 

p c 

PC 

composite properties array 

P/ 

PL 

ply properties array 

P ^p 

PROP 

string PROP; ply proper- 
ties identifiers 
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Engineering 

symbol 

FORTRAN code 

Comment 

Qf,i,p,r,s 

QF,I,P,R,S 

indices to print out 
string PROP 

R 

R 

ply transformation matrix 

*Vi 

RS 

ply transformation matrix 

S j 

SS 

relative rotation matrix 
(between plys) 

S /11T etc - 

P L( 5 1 to 59,1) 

ply limit stresses 

V 

TL 

ply thickness 

VC 

VC 

factor representing void 
effect on strain magni- 
fication factor 

w cb 

WXX 

composite local curvatures 
relative to the structural 
axes 


x c CTL ply coefficient of thermal 

expansion three- 
dimensional in PL(24 to 
26) 

xf ,, m VAF,VAL,VAM filament, ply, and matrix 

’ ’ thermal coefficient of 

expansion 

/? e , y G VCF correlation factors for 

ply thermoelastic proper- 
ties and strain magni- 
fication factors 
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Engineering FORTRAN code 

symbol 


Comment 


A h bta 

H s BET 


correlation factors for 
ply heat conductivities 


correlation factors for 
ply strength 


6 e 


PL(8,I) 


interply layer thickness 


£ csx 


UX 


reference plane membrane 
strains 


EPS,PL(64 to 66,1) ply strains 



THCS 


angle between composite 
material axes and struc- 
tural axes 


i » *Vc 


THLC 


angle between ply material 
axes and composite axes 


v f 12 , m 9 ml2 


NUF12,NUL12,NUM12 filament, ply, and matrix 

Poisson's ratio 


PI 


constant 


RHOF, RHOM, RHOL filament and matrix weight 

density 

SIGL,PL(67 to 69,1) ply stresses 
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